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ABSTRACT 

We estimate the distribution of ejection velocities for the known population of high 
galactic latitude runaway stars. The initial sample is a collection of 174 early-type stars 
selected from the literature. The stars are first classified according to their evolutionary 
status in order to obtain a homogeneous sample of 96 genuine main sequence stars. 
Their present velocities and flight times are then estimated using proper motion data 
from various astrometric catalogues (including Tycho-2, UCAC2, and USNO-B) and 
the ejection velocities are computed by tracing their orbits back in time, based on a 
galactic potential. The potential used is constructed from a mass density model chosen 
to fit the most recent observational constraints. 

We find evidence for two different populations of runaway stars: a "high" veloc- 
ity population, with a maximum ejection velocity of about 400 - 500 kms^"'^, and a 
"low" velocity population, with a maximum ejection velocity of about 300 kms^^. 
We argue that the observed limit of 500 kms~^ and the bimodality of the observed 
ejection velocity distribution are natural consequences of the so-called Binary Ejection 
Mechanism. We discuss the connection between the "high" velocity population and 
the so-called hypervelocity stars, showing how previously studied hypervelocity stars 
are consistent with the results obtained. 

We also find that some stars that were once thought to be best explained as being 
formed in the Halo are compatible with a runaway hypothesis once proper motions 
are included in the analysis. However, three stars in the selected sample appear to 
be inconsistent with ejection from the galactic disc. Possible scenarios are discussed, 
including a possible formation in the galactic halo. 
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1 INTRODUCTION 

Star formation in our galaxy is believed to be confined to the 
star forming regions of the disc. "Runaway" stars are young, 
blue early-type stars that have been observed outside these 
star forming regions (OB associations, open clusters) and 
have usually kinematics different from typical early-type, 
main sequence stars. A surv ey of blue stars performed by 
iGreenstein fc Sargent! (|l974l ) revealed the presence o f these 
main sequence stars in the Halo. In a later work, iTobinI 
argued that a low-resolution spectrum alone does not 
permit us to distinguish normal luminosity main sequence 
stars from the later hot evolutionary stages of low mass stars. 
Other explanations for the presence of runaway stars in the 
Halo were proposed: formation in situ in the Halo, and ejec- 
tion from the disc by a powerful ejection mechanism. Two 
ejection mechanisms have been proposed: 
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(i) the binary ejection mec hanism (BEM): this was 
first proposed by iBlaauwl l|l96ll ) to explain the ejection of 
runaway O and B stars out of the galactic plane. In this sce- 
nario the secondary star of a close binary receives its ejec- 
tion velocity when the primary explodes as a supernova. Be- 
cause one of the stars explodes, this model predicts that run- 
away stars should never be found in binar ies (composed of 
two main sequence stars). Calculations bv lPortegies ZwartI 
(2000) predict a negative mass-ejection velocity correla- 
tion (secondary stars with lower masses receive the greatest 
kicks) and a maximum ejection velocity of < 300 kms~^ . 
Similar result s were found b y Leonard fc DewevI (|l993h . 
More recently, IPrzvbilla et al.l (|2008l ) estimated that veloci- 
ties up to ~ 400 kms~^ are possible in binaries containing 
an early B and a Wolf-Rayet star. 

(ii) the dynamical ejection mechanism (DEM): this 
was first proposed by Poveda, Ruiz & Allen (1967) as an 
alternative to produce runaway stars. Dynamical interac- 
tions between stars inside young, open clusters can give 
large kicks to one or both stars involved in a collision, i.e. a 
close encounter. The large ejection velocities are achieved 
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most efficiently by collisions of two close binaries (since 
they hav e larger cross-sections for the collision). Calcula- 
tions by [Leonard fc Duncan! (|l990l ) predict a binary frac- 
tion among runaway stars of about 10 per cent, a nega- 
tive mass-ejection velocity correlation, and a maxi mum ejec- 
tion v elocity of < 200 kms~^. Later simulations bv lLeonardI 
(|l99ll ). where a much greater number of experiments were 
conducted, revised the value of the maximum ejection ve- 
locity to < 400 kms^^, although very rare events produc- 
ing velocities up to < 1000 kms~^ are possible (by set- 
ting the param eters a posteriori to maximize the velocity). 
More recently, iGvaramadze. Gualandris fc Portegies ZwartI 
(2009) have shown that collisions between binaries and very 
massive stars {M > 50 M0) can also eject stars with veloc- 
ities up to 300 - 400 km s"^ 

Both mechanisms were found to operate in natur e in 
a study by iHoogerwerf. de Bruiine fc de Zeeuwl l|200ll 'l by 
tracing back in time the orbits of runaway stars to their par- 
ent clusters. Two main sequence runaway stars were traced 
back to the same region of the Orion association Ori OBI, 
which is evidence for the DEM. In contrast, a runaway star 
and a pulsar were traced back to the same region of the 
Sco OB2 association, which is evidence for the BEM. This 
study used proper motion (high precision astrometry was 
made available by Hipparcos) and radial velocity data to 
trace back in time the orbits of runaway stars to the par- 
ent cluster. This method makes it possible to estimate the 
moment in time when the star left the cluster and also the 
velocity at that instant (the ejection velocity). Since the or- 
bits are computed in a realistic galactic potential, actual 
estimates of these quantities were obtained instead of lower 
limits as was done in previous studies which used only radial 
velocity information and/or ignored the galactic potential. 
In this work we use a similar method but we apply it to 
distant stars. 

The predicted ejection velocity distribution is similar 
in both models but there are some differences, in particular 
the DEM predicts an enhanced high velocity tail. Studies 
in the past have not been able to properly constrain the 
ejection velocity distribution for two reasons: small, biased 
samples, and lack of proper motion measurements for the 
more distant stars. Since the more distant stars correspond 
to the higher velocities, this second aspect is of particu- 
lar relevance. This is made clear if we remember that the 
more systematic st u dies to date were base d on Hipparcos 
data iMartinl [200^ : I Allen fc KinmanI booj ). with a limit- 
ing magnitude oi V ~ 12.4 and complete up to 1/ ~ 9. 
On the other han d, studies based on t he Palomar-Green 
jCreen et al.ll 19861 ) and Edinburgh-Cape (|Stobie et al.lll997l ') 
surveys reach fainter magnitudes but many stars in these 
studies have unreliable ejection velocity estimates due to 
the lack of proper motion measurements. 

Evidence in recent studies for a link between run- 
away stars and hypervelocity stars has been mounting up. 
Hypervelocity stars are a different class of stars moving 
with even higher velocities than runaway stars, in un- 
bound orbits. They are generally believed to be ejected by 
the S upermassive Black Hole in the centre of the M ilky 
Way (H rown. Geller, Kenvon. Kurtz fc BromlevI [2007bl ) . A 
systematic search for hypervelocity star s under taken by 
iBrown. Geller. Kenvon. Kurtz fc BromlevI l|2007al ) resulted 



in the finding of seven stars still bound to the galaxy. Al- 
though these stars were still interpreted as hypervelocity 
stars (based on their high spatial velocities and large dis- 
tances) and consequently assumed to have been ejected from 
the centre of the galaxy, it is possible that at least some of 
them were ejected far from the galactic centre as well, since 
no dynamical analysis was performed to verify their places 
of origin. The hypervelocity stars HD 271791 and HIP 60350 
were determined to have been ejected from regions far 
from the galacti c centre and with velo cities compatible with 
present models jHeber et al.ll2008| andlPrzvbilla et aLlbOOSl ) 
and for HIP 60350 ijlrrgang et al.|[20m 



The possibility that some young stars that are observed 
in the halo are not consistent with an ejection from the 
galactic plane scenario is still under debate. This is usu- 
ally interpreted to mean that they were formed in situ in 
the halo. This possibility has been investigated in many 
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studies (cf. iHamblv et alll996l and references therein), how- 



ever usually the conclusion is that, once the full kinemati- 
cal information is taken into account, it is not possible to 
rule out the runa way hypothesis . A g ood example is the 
star PHL 346 (cf. |k eenan et all Il986l ') which was shown 
to be consistent wi th ejection from the galactic plane by 
iRamspeck. Heber fc M oehlen (|200ll ). however new candi- 
dates were presented in the same paper, in particular the 
stars PHL 159, PG 1511+467, SB 357, and HS 1914-^7139. 



The purpose of our study is to perform a systematic 
and homogeneous analysis of the sample of known runaway 
stars at high galactic latitudes. This will overcome the ma- 
jor limitations of the previous studies: lack of proper motion 
data in many cases, no error analysis, and limited range 
in apparent m agnitude s /heigh t ab ove the galactic plane. In 
particular, the iMartinI (|200d ') and lAllen fc KinmanI (|2004l 'l 
studies excluded almost all stars in the high velocity tail of 
ejection velocity distribution due to the restriction to bright 
stars, whereas the several s tudies on the Palo mar-Green and 
Edimburgh-Cape surveys (|Stobie et al.lll997l ) lacked proper 
motion data for many stars, although they cover a greater 
range in brightness (and distance). Our aim was to pro- 
duce more solid constraints to the theoretical models by es- 
timating the ejection velocity distribution, shed light on the 
relationship between runaway and hypervelocity stars, and 
investigate the possible existence of candidates for star for- 
mation in the halo. 

This article is organised as follows: in Section [2] we de- 
scribe the selected sample of runaway stars candidates, in 
Section [3] we describe the criteria and method used to clas- 
sify the evolutionary status of the stars in the initial sample, 
including an extensive discussion of the influence of rota- 
tion, ending with a short description of the final sample of 
main sequence stars. In Section |4] we describe the method 
used to compute the orbits of the sample of main sequence 
stars, using the position and full space velocities as input, 
in Section [5] we present our results regarding the ejection 
velocities and flight times, where we try to make a case for 
the existence of an upper limit to the ejection velocity and a 
link between runaway stars and low-velocity "hypervelocity" 
stars, finally we summarize all results in the final Section |6l 
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2 SAMPLE 

Candidate runaway stars were compiled from previous stud- 
ies. We only selected stars for which a spectroscopic analysis 
is published, so we could be sure of their main sequence sta- 
tus and know their atmospheric parameters. The quality of 
the analysis for the different samples is not homogeneous so 
we found it convenient to separate the list of candidate stars 
in two groups. We have included in the first group (Group 
A, Table [l]) the samples selected from the Palomar-Green 
(PG) survey (,Saffer et al... 1997), the Edinburgh-Ca pe (EC) 
survey, and the sample from lRamspeck et al.l 1 200l[ ). In the 
second group (Group B , Tabled we have included the sarn - 
ples from t he pap ers bv lConlon et al.l (|l990l ). lMartinl (|2004 ). 
and lBehJ (|2003bl ). 

The " c omple te" sample of 28 stars selected by 
ISaffer et all (|l997l ) from the PG survey and the sub-samples, 
includin g a "c o mplet e " sample of 13 s tars, selected by 
[RoUcstoi iet al.l (|l997l ). lMagee et al.l (|200ll ). and lLvnn et~aLl 
(2004) from the EC survey, constitute our main sources of 
candidates since they provide good coverage (in a statisti- 
cal sense) of both he mispheres. The PG sample was studied 
in a s eries of papers: 'ifamblv et al.l ()l997l ). iRolleston et al.l 
j 19991 ). and Lvnn ct al. (2004). These studies performed 
high resolution spectroscopic analyses, including the deter- 
mination of radial velocities and atmospheric parameters, 
and detailed abundance analyses. T he EC sample w a s sim i- 
larly studied i n a se ries o f papers by IRolleston et al] ([l993) , 
iMagee et al] (|200ll ) and iLvnn et al.l (|2004l '). Wehave also 
included candidates foun d in fo ur other sources. The sam- 
ple from lRamspeck et al.l (|200ll ) which is of special interest 
because it includes many candidates at high distances (up 
to 7kpc) from the galactic plane (if they are indeed on the 
main sequence). An other important source is the study by 
IConlon et al] l|l992D which contains seven candidates, many 
of which appear to be at very high distances from the plane 
(more than 4kpc). 

The sample of IConlon et al.l l|l990l ) provides 32 can- 
didates which are also part of the Hipparcos cata- 



logue. These 32 stars were studied in a series of papers 
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: Keenan & Duftonll 19831: 


iKeenan, Brown & Lennon 1986;: 


Conlon et al.l 1 19881). Fi- 



Hooi) and iBehJ (|2003bh . iMartinI (|2004l ') performed an ex- 
tensive analy sis of i ts cand idates. It should be noted that 
the study by iBehJ (|2003bD focuses on the distribution of 
the projected rotational velocity of Blue Horizontal Branch 
stars, hence its atmospheric parameters estimates are not 
appropriate for the parameter range of runaway stars. It is 
important to note that there are overlaps between these dif- 
ferent samples. The total number of initial candidate stars 
was then 174. 



3 EVOLUTIONARY STATUS 

In principle, it should be possible to identify main sequence 
stars from the position they occupy on a Tcff — logy di- 
agram, however the region of the diagram which corre- 
sponds to O and B stars is crossed by low mass st ars in 
post main sequence evolutionary stages l|Tobin|[l987l ). Low 
mass stars (initial mass M < 2 Mq) will evolve to the 



horizontal branch after the Helium flash. Stars in the hori- 
zontal branch phase burn Helium in the core and Hydro- 
gen in a shell. The effective temperature of a star when 
it enters the horizontal branch will be directly linked to 
the mass lost during the red giant phase. The Horizontal 
Branch (hereafter HB) stars which occupy the region of 
interest (10000 K < TcB < 30000 K) correspond to the 
hotter end of this mass sequence and are called blue hor- 
izontal branch stars (hereafter BHB stars). The post-HB 
evolution of these stars depends on the mass of the enve- 
lope which determines the strength of th e Hydrogen burning 
shell |d orman. Rood fc 0'Connelllll993D : when the envelope 
mass is greater than a given critical mass, the star will evolve 
to the AGB, and after a period of rapid mass loss, it will 
enter the Post-AGB (hereafter PAGB) phase; on the other 
hand, when the envelope mass is less than the critical value, 
the star will either not reach the tip of the AGB, or not even 
enter the AGB, staying hot until it enters the white dwarf 
cooling sequence. Since BHB stars have envelopes with small 
masses they correspond to this second group. In Fig. [T] we 
have plotted a T^a — log g diagram whore theoretical tracks 
for main sequence (her eafter MS) (S chall er e t al. 1992) and 
low mass PAGB st ars llSchonberneij|l979l ) and the zero age 
horizontal branch (|Dorman et al.lll993l ) are shown, giving 
an indication of the regions occupied with stars in different 
evolutionary stages. 

Stars in our sample were classified as main sequence 
stars or old evolved stars based on their position on the 
Teft — log g diagram, abundance pattern, projected rotation 
velocity, and parallax. In Fig. [TJ we show a — log g dia- 
gram with the stars that were selected from the initial sam- 
ple. These are the stars we believe are on, or near, the main 
sequence and whose orbits were computed, as described in 
Section ID 

3.1 Atmospheric parameters 

Temperatures and gravities were computed from 
Stromgren uvbyj3 photometry l|Moon fc Dworetskvl 
19851 calibration; recalibrated a n d im plemented by 



Napiwotzki. Schonberner fc Wenskd Il993h whenever it 



was available. The main sour ces of photometry were 
iHauck fc MermilliodI (|l998l ) and iMoonev et all (|2000D . the 
latter concentrating only on stars from the PG survey. 
The temperature and gravity estimates obtained in spec- 
troscopic studies were also considered, whenever available. 
These estimates always assume LTE mo del atmospheres 
(with the exception of iLvnn et al. I l2004h . There are two 
groups of estimates: those based on medium-resolution 
spectroscopy and estimates based on high-resolution spec- 
troscopy. Those based on high resolution spectroscopy were 
preferred, followed by the estimates obtained from the 
Stromgren uvhyj3 photometry. The assumed errors on \ogg 
(using the /3 filter) and T^s photometric estimates were of 
0.2 dex and 10 per cent, respectively. The errors on \ogg 
and Tefi estimates taken from the literature range between 
~ 0.1 — 0.25 and ~ 10 — 20 per cent, respectively. iMartinI 
(l2004h estimated \ogg by searching the closest match in a 
model grid. This method is of low precision since it can 
easily produce large errors, potentially twice the size of the 
separation between grid points which was 0.25 dex. Results 
from other investigations were preferred if available. The 
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Table 1. Group A stars, includin g stars from the Palomar-Green (PC) survey l|SafFer et al.lll997^ . the Edinburgh- 
Cape survey, and the sample from lRamsoeck et al.l (2001). The symbol inicl^ no data was available, whereas 
a "(?)" is used to distinguish cases where a MS status is less well established. One / corresponds to 40 kms"-*^ 
V sinj < 70 kms~^, two / correspond to 70 kms~^ ^ v sini < 150 kms~^, and three / to sini > 150 kms~^. 



Name 


Alternative 


Parallax 


u sin i 


Consistent 


Inconsistent 


Verdiet 


References 




name 


< 2 o- 




abundances 


with PAGB nature 






PG 0122 + 214 




7 




Yes 


Yes 


MS 


(6), (7) 


PG 1511 + 367 




7 




Yes 


Yes 


MS 


(6) 


PG 1533 + 467 




7 


^^ ^^ ^^ 


Yes 


Yes 


MS 


(6) 


PG 1610 + 239 




7 




Yes 


Yes 


MS 


(6) 


PHL 159 




7 




7 


Yes 


MS (?) 


(6) 


PHL 346 




7 




Yes 


No 


MS 


(6) 


SB 357 




7 




Yes 


Yes 


MS 


(2), (6) 


BD -15 115 




YCB 




Yes 


Yes 


MS 


(2), (5), (6), (8) 


HS 1914+7139 




7 




Yes 


Yes 


MS 


(6) 


PG 0009 + 036 




7 


■/■/■/ 


7 


Yes 


MS 


(3) 


PG 0855 + 294 




7 


■/ ■/ 


Yes 


Yes 


MS 


(3), (7) 


PG 0914+001 




7 


■/■/■/ 


7 


Yes 


MS 


(3) 


PG 0934+145 




7 




No 


Yes 


Non-MS 


(3) 


PG 0936 + 109 




7 




No 


Yes 


Non-MS 


(3) 


PG 0954+049 




7 


■? 


No 


Yes 


Non-MS 


(3) 


PG 0955 + 291 




7 


■/■/■/ 


Yes 


Yes 


MS 


(3) 


PG 1011 + 293 




7 




No 


Yes 


Non-MS 


(3) 


PG 1205 + 228 




Ygb 


•/ ^ •/ 


Yes 


Yes 


MS 


(1),(3),(7) 


PG 1209 + 263 




7 


■/ ■/ 


No 


Yes 


MS 


(3) 


PG 1212 + 369 




? 


? 


No 


? 


Non-MS 


(3) 


PG 1213 + 456 




? 


? 


No 


Yes 


Non-MS 


(3) 


PG 1243 + 275 




? 


7 


No 


No 


Non-MS 


(3) 


PG 1310 + 316 


PB 3408 


? 


? 


No 


Yes 


Non-MS 


(3) 


PG 1332 + 137 


Feige 84 


Yes 


vv 


Yes 


Yes 


MS 


(3), (7), (8) 


PG 1351 + 393 


PB 890 


? 


7 


No 


Yes 


Non-MS 


(3) 


PG 2111 + 023 




? 




Yes 


Yes 


MS 


(3) 


PG 2120 + 062 




? 




Yes 


No 


Non-MS 


(3) 


PG 2128 + 146 




? 




No 


No 


Non-MS 


(3) 


PG 2134+049 




? 


_ 


No 


No 


Non-MS 


(3) 


PG 2146 + 087 




? 


_ 


No 


Yes 


Non-MS 


(3) 


PG 2159 + 051 




? 




No 


Yes 


Non-MS 


(3) 


PG 2214+184 




7 




No 


Yes 


Non-MS 


(3) 


PG 2219 + 094 




7 


^ ^ ^ 


Yes 


Yes 


MS 


(3), (6), (7) 


PG 2229 + 099 




7 




No 


Yes 


MS (?) 


(3) 


PG 2237+178 




? 


? 


7 


? 


7 


(3) 


PG 2345 + 241 




? 




Yes 


Yes 


MS 


(3), (7) 


PG 2351 + 198 




7 


_ 


No 


Yes 


Non-MS 


(3) 


PG 2356 + 167 




7 




No 


Yes 


Non-MS 


(3) 


EC 04420-1908 




7 




Yes 


Yes 


MS 


(5) 


EC 01483-6804 


CPD -68 91 


7 




No 


No 


Non-MS 


(5) 


EC 05515-6231 


HD 40031 


7 


■/ 


No 


Yes 


MS (?) 


(5) 


EC 06012-7810 




? 




No 


Yes 


Non-MS 


(5) 


EC 09470-1433 




7 




Yes 


Yes 


Non-MS 


(5) 


EC 19071-7643 


HIP 94513 


Yes 




Yes 


Yes 


MS 


(5) 


EC 19337-6743 


HD 184308 


No 


■/■/■/ 


Yes 


Yes 


MS 


(5) 


EC 19476-4109 


HD 187311 


YcB 


■/ ■/ 


Yes 


Yes 


MS 


C5) 


EC 19489-5641 


BPS CS 22896 -0165 


7 




7 


? 


7 


(5) 


EC 19490-7708 




7 




No 


Yes 


Non-MS 


(5) 


EC 19579-4259 




7 




No 


Yes 


Non-MS 


(5) 


EC 19586-3823 


CPD -38 7924 


7 


-/ -/ -r 


Yes 


Yes 


MS 


(5) 


EC 19596-5356 




7 


■/■/■/ 


Yes 


Yes 


MS 


C9) 


EC 20011-5005 




7 




Yes 


Yes 


MS 


C5) 


EC 20089-5659 


HD 191466 


7 




Yes 


Yes 


MS 


(5) 


EC 20104-2944 




? 




Yes 


Yes 


MS 


(5) 


EC 20252-3137 




7 


•/ 


No 


Yes 


MS (?) 


(5) 


EC 20485-2420 




7 




No 


No 


Non-MS 


(5) 


EC 03240-6229 




7 


•/ •/ >r 


Yes 


Yes 


MS 


(4) 


EC 03462-5813 




Yes 


■/■/■/ 


Yes 


Yes 


MS 


(4) 


EC 05229-6058 




7 




No 


No 


Non-MS 


(4) 


EC 05438-4741 




7 




Yes 


Yes 


MS 


(4) 


EC 05490-4510 




7 




No 


Yes 


MS (?) 


(4) 


EC 05515-6107 




7 




Yes 


Yes 


MS 


(4) 


EC 05582-5816 




7 


■/■/■/ 


Yes 


Yes 


MS 


(4) 


EC 06387-8045 




7 




Yes 


Yes 


MS 


(4) 


EC 09414-1325 




7 


■/■/■/ 


Yes 


Yes 


MS 


(4) 


EC 09452-1403 




7 




7 


Yes 


MS (?) 


(4) 


EC 10087-1411 




7 


■/■/■/ 


Yes 


Yes 


MS 


(4) 


EC 10500-1358 




7 


•/ •/ 


Yes 


Yes 


MS 


(4) 


EC 10549-2953 




7 




Yes 


Yes 


MS 


(4) 


EC 11074-2912 




7 




No 


Yes 


Non-MS 


(4) 


EC 13139-1851 




7 




Yes 


Yes 


MS 


(4) 


EC 20140-6935 




Yes 


>r 


Yes 


Yes 


MS 


(4) 


EC 20153-6731 




7 




Yea 


Yes 


MS 


(4) 


EC 20292-2414 




Yes 




Yes 


Yes 


MS 


(4) 


EC 20411-2704 




7 




Yes 


Yes 


Non-MS 


(4) 


EC 23169-2235 




7 




Yes 


Yea 


MS 


(4) 



Mam Jeck el al.| BMH) ; 'i W^ liiuijU ); (Sl IMartij liwa) ; ( 9 1 ILynn et-ST fajiM . 



values for the stars classified as main sequence are shown in 
Appendix [Bl 

The position of a star in the TeS — logp diagram can 
often be used to rule out a PAGB nature. Indeed, as shown 
in Fig. [ll Tcft and log g of most runaway stars are only con- 
sistent with either a MS or BHB evolutionary status. In 
Fig. [T] we can see that the BHB stars populate only a nar- 
row strip in the Teff — logp diagram, however this phase 
is relatively long lived (when compared with the lifetime 
of OB main sequence stars), lasting about 100 Myr, which 



means this is the region of the Tch — log g diagram where we 
expect most contamination by low mass evolved stars. How- 
ever, most stars classified as MS are in the region between 
the PAGB tracks and the terminal BHB. Given that BHB 
stars spend only ~ 10 Myr in the post-BHB phase (which is 
sometimes called "AGB-mangue" , as these stars do not as- 
cend the AGB), the contamination in this region is expected 
to be minimal. The number of selected stars in the region 
where the PAGB tracks intersect the main sequence is even 
smaller. As an aside, we note that the lack of stars in the 
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Table 2. Group B stars, including the samples from the papers of IConlon et al ] l|l99d ). lMartinl l|2004l 'l. and iBehJ 
ll2003bl'l . The symbols have the same meaning as in Table [T] 



Name 


Alternative 
name 


Parallax 

< 2 


V sin i 


Consistent 
abundances 


Inconsistent 
with PAGE nature 


Vcrdiet 


Refcrenecs 


PB 5418 




? 




Yes 


Yes 


MS 


(2) 


Ton S 195 


SB 463 


? 


- 


No 


Yes 


Non-MS 


(2) 


Ton S 308 




? 




Yes 


Yes 


MS 


(2) 


PHL 2018 




? 


■/■/■/ 


Yes 


Yes 


MS 


(2) 


BD -2 3766 


HD 121968 


Yes 




Yes 


Yes 


MS 


(2), (8) 


BD +00 0145 




7 


- 


No 




Non-MS 


(7) 


BD +36 2242 




Yes 




7 


Yes 


MS (?) 


(7) 


HD 7374 




No 


- 


7 


Yes 


Non-MS 


h) 


HD 27295 




No 


- 


7 


Yes 


Non-MS 


(7) 


HD 128801 




No 


- 


No 


Yes 


Non-MS 


h) 


HD 135485 




No 


- 


No 


Yes 


Non-MS 


(7) 


PG 1530 + 212 




? 




7 


Yes 


MS (?) 


(7) 


HIP 1241 


HD 1112 


Yes 




No 


Yes 


MS (?) 


(7), (8) 


HIP 1511 




Yes 


- 


No 


Yes 


Non-MS 


(8) 


HIP 1904 


HD 1999 


Yes 


? 


Yes 


Yes 


MS 


(1) 


HIP 2702 


HD 3175 


Yes 


- 


Yes 


Yes 


MS 


(1),(5) 


HIP 3812 


JL 212 


Yes 


^ ^ ^ 


Yes 


Yea 


MS 


(1),(5) 


HIP 6419 


HD 8323 


Yes 


7 


7 


Yes 


7 


(1).(8) 


HIP 11809 


Fcigc 23 


Yes 


^ ^ ^ 


Yes 


Yea 


MS 


(1),(8) 


HIP 11844 


HD 15910 


Yes 




Yes 


Yes 


MS 


(s) 


HIP 12320 


Feige 25 


Yes 


y y y 


Yes 


Yea 


MS 


(1),(8) 


HIP 13800 


Fcigc 29 


Yes 


— 


Yes 


Yea 


MS 


(2), (7) 


HIP 15967 


HD 21305 


Yes 


- 


No 


Yes 


Non-MS 


(8) 


HIP 16130 


HD 21532 


Yes 




Yes 


Yes 


MS 


(1),(8) 


HIP 16466 


HD 21996 


No 


- 


Yes 


Yea 


Non-MS 


(1) 


HIP 16758 


HD 22586 


Yes 




No 


No 


MS (?) 


(1),(5) 


HIP 28132 


HD 40267 


Yes 




Yes 


Yea 


MS 


(8) 


HIP 37903 


BD +61 996 


Yes 




Yes 


No 


MS 


(8) 


HIP 41979 




Yes 


- 


No 


Yea 


Non-MS 


(8) 


HIP 45904 


HD 233622 


Yes 


-/-/-/ 


No 


Yea 


MS (?) 


(7), (8) 


HIP 48394 


HD 237844 


Yes 




No 


Yea 


MS (?) 


(8) 


HIP 50750 


BD +16 2114 


Yes 


7 


7 


Yea 


7 


(8) 


HIP 51624 


HD 91316 


Yes 


? 


7 


No 


7 


(1) 


HIP 52906 


BD +38 2182 


Yes 




Yes 


Yea 


MS 


(1),(8) 


HIP 55051 


HD 97991 


Yes 




? 


No 


MS (?) 


(1) 


HIP 55461 


Fcigc 40 


Yes 




Yes 


Yea 


MS 


(1),(7),(8) 


HIP 56322 


HD 100340 


Yes 




Yes 


Yea 


MS 


(7), (8) 


HIP 58046 


HD 103376 


Yes 


y y y 


Yes 


Yea 


MS 


(1)>(7),(8) 


HIP 59067 


HD 105183 


Yes 




Yes 


Yea 


MS 


(1)>(7),(8) 


HIP 59955 


HD 106929 


Yes 




Yes 


Yea 


MS 


(1),(8) 


HIP 60578 


BD +49 2137 


Yes 


- 


7 


Yea 


7 


(7), (8) 


HIP 60615 


BD +36 2268 


Yes 




No 


Yea 


Non-MS 


(1)>(7),(8) 


HIP 61800 


HD 110166 


Yes 


y y y 


Yes 


Yea 


MS 


(8) 


HIP 65388 


PB 166 


Yes 


- 


No 


Yea 


Non-MS 


(8) 


HIP 69247 


HD 123884 


Yes 


- 


No 


Yea 


Non-MS 


(8) 


HIP 70275 


HD 125924 


Yes 




7 


Yea 


MS (?) 


(1)>(7),(8) 


HIP 71667 


BD +20 3004 


Yes 






Yea 


MS (?) 


(7), (8) 


HIP 75577 


HD 137569 


Yes 




Yes 


Yea 


Non-MS 


(7), (8) 


HIP 76161 


HD 138503 


Yes 




Yes 


No 


MS (?) 


(8) 


HIP 77131 


HD 140543 


Yes 


y y y 


No 


No 


MS (?) 


(8) 


HIP 77716 


BD +33 2642 


Yes 




No 


No 


Non-MS 


(7), (8) 


HIP 79649 


HD 146813 


No 




Yes 


Yea 


MS (?) 


(1),(8) 


HIP 81153 


HD 149363 


Yes 




No 


No 


MS (?) 


(8) 


HIP 82236 


BD +13 3224 


Yes 


7 


7 


Yea 


7 


(8) 


HIP 96130 


HD 183899 


No 




Yes 


Yea 


MS (?) 


(8) 


HIP 98136 


HD 188618 


Yes 




Yes 


No 


MS (?) 


(8) 


HIP 104931 




Yes 


7 


? 




7 


(8) 


HIP 105912 


HD 204076 


Yes 




Yes 


No 


MS (?) 


(1) 


HIP 107027 


HD 206144 


Yes 




Yes 


No 


MS (?) 


(1),(8) 


HIP 108215 


HD 208213 


Yes 


y y y 


Yes 


Yea 


MS 


(1) 


HIP 109051 


HD 209684 


Yes 




Yes 


Yea 


MS 


(8) 


HIP 111396 


HD 213781 


Yes 




Yes 


Yea 


MS 


(1),(7),(8) 


HIP 111563 


HD 214080 


Yes 




Yes 


No 


MS (?) 


(1) 


HIP 112790 


HD 216135 


Yes 


■/ ■/ 


Yes 


Yea 


MS 


(1),(8) 


HIP 113735 


HD 217505 


Yes 


7 


Yes 


Yea 


MS 


(1) 


HIP 114569 


HD 218970 


Yes 




No 


Yea 


MS (?) 


(1),(8) 


HIP 114690 


HD 219188 


Yes 




Yes 


Yea 


MS 


(1) 


HIP 115347 


HD 220172 


Yes 




Yes 


No 


MS (?) 


(1) 


HIP 115729 


HD 220787 


Yes 




Yes 


Yea 


MS 


(1)>(7),(8) 



(i)ES3S 



et all ^199qi: r2^ |Conlon et all 119931; f31 | Satter et all flE^Tjl: f4 ^ IRollestc 
IRamBpeck et all 120011 1; (TT TBehrl I2003bl l; (SI IMartiiJ I2004II ; l91 ILvn 



] {1991); 

r i200'yi . 



(Sl IMagee et all t200Jl ; (6) 



PAGE region is due to the small number of O stars in our 
sample. In summary, BHB and post-BHB stars occupy the 
region of the Tetr — logp diagram between the ZAHB and 
the PAGB tracks, however they spend 90 per cent of their 
lifes between the ZAHB and the TAHB, and only 10 per 
cent between the TAHB and the PAGB track. Hence the 
contamination after the TAHB should be minimal. 

3.2 Abundances 

An abundance analysis was, in general, part of all the spec- 
troscopic studies. This analysis uses the respective atmo- 
spheric parameters estimates and assumes an LTE atmo- 
sphere. The abundance pattern permitted us to distinguish 
MS stars from BHB stars, as the atmospheres of the lat- 
ter are dominated by diffusion and show strong deviations 



from the approximately solar mix seen in MS stars. In par- 
ticular, there is usually an He depletion and depletion or 
enhancement of metals, depending on the evolutionary his- 
tory, and balance between the effects of diffusion and levi- 
tation of the heavier metals due to radiation pressure (|Behj 
l2005h . In the case of PAGB stars, some elements can be 
modified by dredge ups, but the absence of strong modifi- 
cation of the abundance pattern is usually not conclusive. 
The abundances found in the literature were c ompared with 
the n ormal abundances for B stars, found in iKilianI (Il992l . 
11994 ). A star with abundances that were not consistent with 
normal abundances, within errors, was marked as such. 
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reff(K) 

Figure 1. T^s — logg diagram. The white squares with error bars represent stars classified as main sequence. The black diamonds 
repre sent stars classified as evolved stars, mostly blue horizontal branch. Also shown are theoretical tracks for masses in the range 3-20 
Mq l lSchaller et al.lll992l) as solid lines. The zero age horizontal branch (ZAHB) is shown as a dashed-dotted line and th e terminal age 
horiz ontal branch (TAHB) is shown as a dashed line, both for a Helium mass fraction of Y = 0.247 and [Fe/H] = -1.48 ||d orman et al.l 
Il993l) . Note that the tracks for the runaway stars are for solar metallicity, in contrast with the (low) metallicity assumed for the HB 
stars, which is appropriate for halo metal-poor stars. The theoretical track of a PAGB star of 0.546 Mq llSchonberner|[l979l) is shown as 
a dotted line. 



3.3 Projected rotation velocity 

A high projected rotation velocity, or vsini, where i is the 
angle between the line of sight and the rotation axis, is an 
excellent indicator of a young age as old evolved stars do not 
rotate as fast as young objects. The maximum rotation ve- 
locity for_ablue HB star appears to be around 30 — 40 kms^^ 
l|Behij|2003bl ). As such, we have considered a star with a pro- 
jected rotation velocity in excess of 40 kms~^ to be likely 
a young object, a star with a projected rotation velocity 
in excess of 70 kms~^ to be very likely a young object, 
and a star with a projected rotation velocity in excess of 
150 kms~^ to be extremely likely a young object. This is 
noted in Tables [1] and [2] by 1, 2 or 3 ticks respectively. 

It is known that a high rotation velocity will make a 
star appear cooler and more luminous (lower surface grav- 
ity impl ying a larger radius ) if observed at high incl i natio n 
angles ( Fremat et al ] l2005l : IWenske fc Schonbernerl Il993t ). 
This effect is sometimes called gravitational darkening and is 
caused by the increasing effective gravitational acceleration 
as a function of the latitude, with the equatorial regions be- 
coming cooler than the polar regions of the s tar. The selected 
samp le contains three Be stars: SB 357 (iRamspeck et al] 
I2OOII). PG 0914-K) 01 (|Rolleston eral]|l999l ) and HIP 3812 
i Magee et aLlliooH ). Typical rotation velocities of a Be star 
are very hig h (on average 88 per cent of the break-up velocity 
according to lFremat et al.l2005l ). we would expect significant 
effects due to gravitational darkening. Indeed, if the mea- 



sured surface gravity for PG 0914+001 is taken at face value 
then the estimated critical velocity {vcrit ~ 273 kms~^) 
will be lower than the measured projected rotation veloc- 
ity (usini = 325 kms~^), which strongly suggests that the 
surface gravity has been underestimated implying a lumi- 
nosity overestimate. This would also explain the very large 
distance derived for the star PG 0914-1-001 of ~ 35 kpc, 
based on a low measured surface gravity. Although the grav- 
itational darkening effect is more important in Be stars, be- 
cause of their high rotation velocity, it should not be ignored 
in nor mal B stars as pointed out by Wenske fc Schonberneij 
(|l993l ). According to lAbt et all l|2002l ) B stars rotate with 
velocities of 40 - 50 per cent of the critical velocity on 
average (depending on the exact spectral type). On the 
other hand, the effect only becomes significant for veloc- 
ities greater than ~ 60 per cent of the critical velocity 
l|Fremat et al.ll2005h . 

We have corrected the Be stars observed tempera- 
tures and sur f ace g ravities using the theoretical models by 
iFremat et all l|2005l ) assuming a rotation velocity of 99 per 
cent of the break-up velocity. Although there is some indi- 
cation that many Be stars are actually rotating slower than 
this (88 per cent i s the most probable value according to 
iFremat et"al]|2005l. and 70 per cent the lower limit accord- 
ing to lEkstrom et al.ll20()8l ) , we know that at least in the case 
of PG 0914-1-001 the rotation must be very close to critical 
because of the high vsini. The case for rotation close to the 
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critical velocity is weaker for the two other Be stars, but 
the same correction is applied due to homogeneity consid- 
erations. Nevertheless, the difference between a correction 
assuming a rotation of 99 and 80 per cent of critical velocity 
is minimal for these two other Be stars (~ 0.1 dex for logg 
and 1000 K for Teff). 

In the case of the normal B stars, the gravitational dark- 
ening effect only becomes significant for velocities greater 
than 60 per cent, as was mentioned previously. Neverthe- 
less, the model for 80 per cent of break-up velocity was pre- 
ferred because some stars have a projected rotation velocity 
value incompatible with a rotation velocity lower than 60 
per cent of the critical velocity. Both approaches produce 
similar results (the corrections for \ogg and Tofi are, on av- 
erage, 0.1 dex and 400 K for the 60 per cent case, and 0.1 dex 
and 485 K for the 80 per cent case). 

In principle we would like to correct the gravitational 
darkening effect for all the stars in the sample, but, since 
the true rotation is unknown for any given star, the cor- 
rection is done in a statistical sense only. The chosen model 
implies a very large true rotation velocity (80 per cent of the 
critical rotation velocity), hence we would overestimate the 
correction for most stars with small v sin i values. To account 
for this problem we have done the correction only for stars 
with V sin i larger than 35 per cent of the break-up velocity. 
Moreover, the correction was not applied in cases where ro- 
tation velocity measurements were not available (HIP 1904, 
HIP 113735 and HD 138503), and when the effective tem- 
perature was much higher than the maximum temperature 
(27000 K) available in the theoretical model grid (HD 140543 
and HD 149363). 



3.4 Selection results 

By applying the selection criteria, we have classified 96 stars 
of the initial sample of 174 as being likely or very likely on 
the main sequence. The remaining 78 stars are most likely 
halo population, evolved stars, mostly stars on the horizon- 
tal branch judged from their position on the T^fi — logg 
diagram. The 96 selected stars have been further classified 
according to how strong their case for being on the main 
sequence is. In the more convincing cases there is good evi- 
dence for a normal abundance pattern, and/or for high rota- 
tion velocity, whereas the less convincing cases either present 
weaker evidence for normal abundance coupled with low or 
non-measured rotation velocity. These two cases correspond 
to the verdicts of MS and MS (?) in Tables [T] and [2] respec- 
tively. All relevant data for this sample of main sequence 
stars is shown in Appendix [Bl 

The selected sample covers a range in magnitudes of 

6.5 < V < 14.5, which corresponds to a range in heights 
above the galactic plane of 0.3 — 30.5 kpc. It is this sample of 
main sequence stars that is the focus of interest throughout 
the remainder of this paper. 



4 ORBITAL ANALYSIS 

The method used to estimate the ejection velocities consists 
in tracing the orbits of the stars backwards in time until the 
first intersection with the galactic plane. The present space 



velocity is obtained from the radial velocity and proper mo- 
tion (given the distance), which are available for every star 
in our sample. Given the velocity and the position of the 
star, we integrated the orbit in the galactic gravitational 
potential after inverting the velocity direction (to go back 
in time). This is esse n tially the m ethod used to good effec t 
by iHoQgerwerf et all (|200ll '). and iRamspeck et al.l l|200j ). 
among others. The final output is then the instant of inter- 
section with the galactic plane - which is equal to the time 
the star spent on the orbit after ejection (or fiight time) 
- and the velocity at that instant (the ejection velocity). 
The orbi ts were computed using the pro gram ORBItO devel- 
oped bv lOdenkirchen fc Broschd l| 19921 '). This program m- 
teg rates the orbits of t est pa rticles in a modified version of 
the lAllen fc SantillanI (|l99ll ) potential, where the disc has 
a scale- length of 3 kpc. The changes made to the original 
galactic potential are detailed in Appendix|X] The program 
takes as input the spatial coordinates and full space velocity 
(in a right-handed galactocentric frame of reference where 
the X axis passes through the position of the Sun, pointing 
in the opposite direction, and Z points to the north galac- 
tic pole). The steps involved in computing these values, in 
chronological order, were: 

• the mass was obtained from the atmospheric parame- 
ters, by interpolating the theoretical tracks shown in Fig. [T] 

• the distance was computed from the visual magnitude 
and the absolute magnitude, which in turn was obtained 
from the mass and atmospheric parameters; 

• the coordinates and velocities in the galactocentric 
Cartesian frame of reference were obtained from the mea- 
sured equatorial coordinates, radial velocities and proper 
motions, and distances. 

Each of these steps was implemented as a separate subrou- 
tine in a FORTRAN program. Errors were propagated using 
a Monte Carlo scheme by assuming a Gaussian distribu- 
tion and choosing a random point from the input parameter 
space on each iteration. All steps will be expanded on in 
later subsections. 

The orbit for each star was integrated 10000 times 
(in the Monte Carlo scheme previously described) using a 
timestep of 0.2 Myr, and the program stopped, on any given 
iteration, when the orbit reached the galactic plane (which 
corresponds to Z = in the galactocentric Cartesian refer- 
ence frame) or when the flight time exceeded 250 Myr (to 
avoid the program running forever in those situations where 
the simulated star did not reach the plane), whichever hap- 
pened first. The velocity of the star when it reached the 
galactic plane, after being corrected for the rotation of the 
galactic disc, was taken as the ejection velocity. The time 
since the beginning of the integration until the crossing of 
the plane was taken as flight time. Instances with the star 
not reaching the galactic plane within 250 Myr were not in- 
cluded in the calculation of velocities and flight times. 

4.1 Masses, ages and distances 

The masses and ages were obtained by int erpolating be- 
tween the theoretical evolutionary tracks of ISchaller et al.l 
(|l992l ) (Schaller tracks). We note that for this reason, 
the age of a star is more properly termed evolution- 
ary age. The tracks were converted to the Tcff — logg 
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plane. A metallicity corresponding to Z = 0.02 (close 
to solar) was assumed. As a test of the uncertainty 
caused by different input physics, the evolutionary ages 
were compared to diff erent estimates obtained using the 
theore tical tracks by ISressan. Fagotto. Bertelli fc Chiosil 
(|l993t ) ( Padova tracks), which are the same as the 
tracks bv lGirardi. Bressan. Bertelli fc Chiosil 12000) during 
the main sequence phase, and to the theoretic al tracks 
by iPietrinferni. Cassisi. Salaris fc Casteiiil (|2004l ') (BaSTI 
tracks). This latter tracks are only available for stars with 
less than 10 M© (note that this covers 80 per cent of the 
sample). In the case of the Padova tracks, the relative dif- 
ference between the two determinations is typically less than 
10 per cent for stars older than 50 Myr and less than 20 per 
cent for stars younger than 50 Myr. Note that the increased 
discrepancy for younger stars can be partly attributed to 
the fact that the two models start at different evolutionary 
stages. In the case of the BaSTI tracks, the relative differ- 
ences in evolutionary age are typically less than 10 per cent 
for all stars (masses lower than 10 M©). In both compar- 
isons, the differences can at least partly be explained by the 
different treatment of the overshoot into the convective lay- 
ers. Note that these differences between sets of theoretical 
tracks are in most cases much smaller than the errors in the 
determination of the evolutionary age, resultant from the 
uncertainty in the effective temperature and surface gravity 
determinations. 



The distance was obtained from the distance modu- 
lus. Absolute magnitudes were computed for given logg, 
Tcft and mass using V band fluxes from ATLAS model 
atmospheres. We note here that the apparent magni- 
tudes were corrected for interst ellar reddening using eit her 
Stromgren uvby(3 photomet ry (INapiwotzki et al.l 19931 re- 
calibration for B stars of the lMoon fc Dworetskvlll985l orig- 
inal calibration) when available, or the reddening maps of 
ISchlegel. Finkbeiner fc Davi s ( 1998). These reddening maps 
include the total galactic reddening in any given direction. 
However we are dealing with stars in the Halo only and most 
reddening occurs in the disc. This was verified by comparing 
the Stromgren uvbyP photometry reddening estimates with 
the Schlegel values, which were in good agreement. 



The spectroscopic distance determination can be tested, 
if accurate trigonometric parallaxes exist. However, even the 
accuracy achievable with Hipparcos is not sufficient for the 
sample of runaway stars. Only seven stars have parallax er- 
rors below 50 per cent and the best case corresponds to an 
error > 30 per cent. However, we could carry out a test of the 
spectroscopic metho d using the sample of ea rly type stars 
from the Table 4 of INapiwotzki et al.l (|l993l ). Parameters 
were derived in a fashion similar to many investigations of 
runaway stars: Tefi was determined from Stromgren photom- 
etry and log g from Balmer line fitting. All these stars have 
accurate Hipparcos parallaxes. We found that the ratio of 
distances determined using our method and the trigonomet- 
ric parallax determination is on average 1.03, with a stan- 
dard deviation of 0.11. This dispersion is of the same order 
of the expected error in distance corresponding to the best 
determinations of logg, which is ~ 10 per cent for an error 
of 0.1 dex. 



4.2 Position and velocity 

Current positions in equatorial coordinates were taken from 
the Hipparcos and the UCAC 2 catalogues, using the 
equinox 2000.0 transformation as given by SIMBAD, and 
then converted to three-dimensional galactocentric Carte- 
sian coordinates, {X,Y,Z). Velocities were computed from 
measured proper motions and radial velocity (relative to the 
local standard of rest - LSR) and transformed to galacto- 
centric Cartesian coordinates by accounting for the motion 
of the Sun and of the LSR. The adopted Sun - Galactic 
Centre distance was 8.0 kpc and the Sun velocity in galac- 
tocentric Cartesian coor dinates was (11.0,5.3,7.0) kms~^ 
iDehnen fc Binnevlirggi ). 

The sources used for the pr oper motions were: the 
UCAC2 catalogue (^charias et al.. 2004 ). the USNO-B cat - 
alogue jMonet et al.ll2003l). theTy cho-2 (|H0g et all [200ol ) 



and Hipparcos (|van LeeuwenllioOT ') catalogue s, the Super 
COS MOS science archive, the SP M catalogue (ICirard et al.l 
l200i), and the NPM2 catalogue (|Hanson et all 120031 ). The 
different proper motion estimates were combined by a 
weighted average following the prescription in Section 2 of 
iPauli et al.l (|2006l ). The values obtained are shown in Ap- 
pendix IB] 

Many proper motion measurements in the USNO-B cat- 
alogue lack estimates of the associated error, or in other 
cases an error of zero is given. To account for these two situ- 
ations we have estimated the typical proper motion error for 
the stars in the sample. This was accomplished by compar- 
ing the USNO-B proper motions with the other catalogues, 
in particular with Tycho-2 and UCAC2, and computing the 
USNO-B error which minimized the weighted (with the er- 
rors) average difference between the proper motions from 
USNO-B and the other catalogues. 

It was determined that the quality of the Super- 
COSMOS proper motion measurements tended to decrease 
with increasing brightness. For this reason, SuperCOSMOS 
proper motions were considered only for stars fainter than 
V — 14. This limit was determined from comparisons with 
all other proper motion sources. 



4.3 Error propagation 

The orbit obtained at the end is a very complicated function 
of the various input parameters, but still it is important to 
have an estimate of the errors in the ejection velocities and 
flight times. An error analysis was ignored in most previous 
studies wit h the notable exception of the studies bv lMartinI 
lloO^) and lRamspeck et al.l 1200 1'). Since it is impossible to 
estimate these errors using the standard error propagation 
formula we implemented a Monte Carlo simulation where 
each input parameter is a random point taken from a Gaus- 
sian distribution which reflects the error distribution of the 
measurement. A further advantage of this procedure is that 
it permits a study of the correlations between the different 
variables. 



5 RESULTS AND DISCUSSION 

The flnal results of the analysis described in Section |4l per- 
formed on the 96 main sequence stars selected from the ini- 
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tial sample of 174, are the flight time (i.e. the time spent in 
the orbit since ejection until the present), the ejection veloc- 
ity (i.e. the orbital velocity for Z = 0, relative to the star's 
LSR at the moment of ejection) and spatial {Xq, Yq) coordi- 
nates of the ejection point, and their associated errors. The 
distributions of these output variables are highly skewed, 
hence the median is used as a measure of central tendency 
value instead of the average, and the percentiles (confidence 
interval bounds) corresponding to 15.9 and 84.1 per cent of 
the distribution (corresponding to ilcr for a Gaussian dis- 
tribution) were used instead of the standard deviation as a 
measure of the error range. 

The model of the galactic potential is an obvious source 
of uncertainty in the final results. A detailed analysis of the 
galactic potential parameter space is outside the scope of our 
study. However we have tested a number of alternative mod- 
els as an attempt to quantify the effect on the flight times 
and ejection velocities, by choosing different, reasonable, pa- 
rameters, starting from the base model (cf. Appendixl^. We 
have tested three modified models: one with a disc scale- 
length shorter by 50 per cent, one with a disc scale-length 
longer by 50 per cent, and one with a total mass higher by 
30 per cent. In all cases the disc component mass, and the 
halo component mass and power law index were changed to 
fit the galaxy's circular rotation curve to the observed value 
of 220kms~^ at the position of the Sun. The plot of the 
relative difference of each model to the default model, as 
a function of height above the galactic plane, is shown on 
Fig. [2l for the flight times, and Fig. [S] for the ejection veloc- 
ities. It is noticeable that the difference between the models 
drops with distance to the disc. It is not surprising that the 
influence of the disc component becomes less significant for 
large distances. Note in particular that the model with a 
heavier halo only affects the outcome because the choice of 
a heavier halo forces a much lighter disc in order to keep 
the local circular velocity unchanged. However, in the outer 
regions, where the disc component loses its influence, the 
results are almost unchanged which suggests that the mass 
of the halo is not crucial for the obtained results, as long 
as the local orbital velocity is kept at ~ 220 km s^^. The 
fact that the results for the stars which are further away are 
mostly insensitive to the choice of potential is particularly 
important, since these are the stars that have the highest 
ejection velocities. The difference is more dramatic for the 
flight times which reach a difference of 25 per cent. Never- 
theless, most cases in both the ejection velocities and flight 
times are contained in the interval 10-15 per cent in relative 
difference. 



5.1 Flight times 

If the disc ejection hypothesis is correct it follows that the 
computed flight times must be less than, or equal, to the 
corresponding evolutionary ages (cf. Section [4. 1|) . The flight 
times of all stars in our sample are plotted in Fig. [l] We 
consider them to be consistent within the errors if the flight 
time estimate minus its error, is equal to or lower than the 
evolutionary age estimate plus its error. We find that the 
flight times are consistent with the computed evolutionary 
ages to within the errors for 83 per cent of the sample and 
to within twice the errors for 97 per cent of the sample (cor- 
responding to only three discrepant objects). Moreover, we 
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Figure 2. Relative difference of the ejection velocities between 
the tested galactic potential models. The circles represent the 
short scalelength disc model, the crosses represent the long scale- 
length disc model, and the squares represent the heavier halo 
model. 
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Figure 3. Relative difference of the flight times between the 
tested galactic potential models. The symbols have the same 
meaning as in the ejection velocities plot. 



note that for most stars Tf ~ Ti, as it is ex pected from the- 
ory fo r both the BES and DES mechanisms. iPortegies ZwartI 
([2Qo3) estimates, for the BES, that a late- type B star (mass 
in the interval 3-5 Mq) will spen d more than 75 per cent o f 
its lifetime as a runaway star, and lLeonard fc DuncanI (|199G| 1 
estimate, for the DES, a time for ejection from a cluster of 
the order of 10 million years, corresponding to more than 
75-80 per cent of the lifetime spent as runaway, in the case 
of late-type B star. The fact that we find this agreement with 
theory suggests that the models are adequate or conversely, 
if we assume the models are correct, that the method used 
to compute the fiight times does not have a fatal flaw. In 
particular, we are led to believe that the assumed galactic 
potential (cf. Section|4]and AppendixIXJ is realistic enough 
for our purposes. This conclusion is supported by the dis- 
cussion in the beginning of the section on the sensivity of 
the orbits to changes in the galactic potential. In spite of 
this overall good agreement of observations and theoreti- 



10 M. D. V. Silva and R. Napiwotzki 




Age (Myr) 

Figure 4. Flight time versus evolutionary ages including error 
bars. 
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Figure 5. Flight time versus evolutionary ages including error 
bars. Only stars deviating more than twice the errors in both 
directions are plotted. 



cal expectation, we find three stars in our sample for wliicli 
tfie computed fiiglit time is in serious disagreement (more 
tfian twice tfie errors in botli directions) witti our estimate 
of tlieir evolutionary age (see Fig. [5} . Taken at face value 
tfiis could indicate that thes e stars were born in the galactic 
halo, a scenario proposed bv lHamblv. Wood. Keenan et aL] 
for the star PHL 346. If confirmed, this could have 
very important consequences for our understanding of star 
formation. However, first we have to consider other possible 
mechanisms and effects capable of extending the lifetime of 
main sequence B stars. These are: 

(i) rotation induced mixing: early-type stars have gen- 
erally high rotation velocities, in particular in the case of 
Be stars which rotate with velocities near the critical ve- 
locity. The centrifugal forces created by the fast rotation 
create extra mixing acting as an extra overshoot diffusing 
the elements produced in the core. This effect would in- 



crease the lifetime of a B star at most by 20-25 per cent 
jMaeder fc Mevnedlioool ). 

(ii) lower than solar metallicity: as described in Sec- 
tion 14.11 ages were computed assuming a metallicity close 
to solar. A significant scatter of metallicities is observed for 
young st ars in the disc. T he metallicity distribution, as de- 
rived bv lFuhrmannI l|2004h from a sample of main sequence 
stars from B to G spectral types within 25 pc of the Sun, 
has a standard deviation of 1.4 times the solar metallicity, 
which is comparable to the radial scatter in the thin disc 
( jCescutti et al.ir2007h . This effect is relevant because stars 
with lower metallicities stay longer on the main sequence, 
and the correspondent ZAMS (zero age main sequence) is 
bluer, implying longer lifetimes for given effective tempera- 
tures and surface gravities, when compared with stars with 
higher metallicity. 

(iii) blue straggler: the DEM may result in the ejection of 
a binary which could then merge after the ejection, creating 
a rejuvenated star analogous to the blue stragglers observed 
in clusters. This situation would be effectively the same as a 
star with a mass equal to the sum of the two stars in the pro- 
genitor binary being ejected from a position higher above the 
galactic plane. Moreover, since the binary stars would have 
a spectral type later than B (to be able to produce a B star 
after the merger) this would increase dramatically the life- 
time of the system. However, the ejection of intact binaries is 
predicted to be a rare e vent and the ejection veloc ities lower 
than about 100 km s~^ (jLeonard fc Duncanlll990l ). More re- 
cently [Peret^ (|2009l ) has demonstrated that this mechanism 
could explain, in principle, the extreme youth of the these 
distant stars, but the problem with the small numbers and 
velocities still remains. 

The alternatives to these mechanisms would be a formation 
in the Halo scenario, or errors in the analysis of the stars' 
atmospheres (which would propagate to the estimated dis- 
tance) and/or errors in the measurement of proper motions. 
The three stars under consideration, that have flight times 
higher than the evolutionary ages within two times the er- 
rors, are: SB 357, EC 20252-3137 and HIP 77131. Ahhough 
a few 2-sigma outliers are expected (about five in a sample 
of 100 objects), these stars have small errors, which makes it 
difficult to explain the discrepancies as statistical fiuke. The 
estimated ages for SB 357, EC 20252-3137 and HIP 77131 
are G.ltH Myr, 12.5lg;4 Myr and 4.5l2:9 ^Y^^ and the 
fiight times 53. 3l^;^ Myr, 45.7l5;^ ^yr and IT.Qtil 
respectively. We have computed the ages the stars would 
have if we take into account the maximum effects of metal- 
licity and rotation. The results are summarized in Table O 
A more detailed discussion of each star follows: 

SB 357 This is a Be star, so its MS status appears to be 
well supported. However a high value for the projected ro- 
tation velocity also means that no metal abundances were 
determined for this s tar, although the Heliu m abundance ap- 
pears to be normal l|Ramspeck et al]|200ll ). The ages given 
in Table |3] are justified by the fact that this star is a fast 
rotator coupled with the fact that there is no evidence for 
Solar metallicity. However, even when both effects are com- 
bined, the difference between ages and fiight times is still 
significant within (the equivalent to) 1.5 a. A different pos- 
sibility is that the gravitational darkening may have been 
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Table 3. The evolutionary ages of SB 357, EC 20252-3137 and 
HIP 77131 if we assume a mass fraction of metals of Z = 0.008, 
and we take into account the effect of rotation. The first col- 
umn gives the ages adjusted for rotation, the second column 
gives the ages for low metallicity, and the third column gives 
the age for metallicity together with effect of rotation. All ages 
are given in Myr. The evolutionary age s for low meta l licites were 
estimated from evolutionary tracks by iFaeotto et al.l l ll99'tf) . We 
have assumed rotation i ncreases the lifetimes by 25 per cent (cf. 
iMaeder fc Mevneti[200ch . 
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overestimated in this particular case. If the correction is ig- 
nored the computed flight time is 65'^l^'^ Myr, whereas the 
age is at most (assuming a low metallicity and adding the 
effect of rotation) 4612 3 Myr, implying that even in this 
situation the difference is still significant within the error 
bounds. We note, however, that since SB 357 is a Be star 
this is a situation where the correction is the most impor- 
tant. As a final test, we compared the age esti mate with the 
one ob tained using the theoretical tracks bv lBressan et all 
l|l993l ), but this increases the discrepancy even more since 
this new determination is a factor of 2 smaller than the orig- 
inal one. 



EC 20252-3137 The classification of this star as being on 
the main sequence appears to be solid. It is rotating with 
a (projected) velocity of 60 kms~^ and has a normal abun- 
dance pattern, although there is some evid ence for slightly 
supersolar metallicities l|Magee et al.ll200ll ). Given this evi- 
dence for high metallicity the assumption of low metallicty 
is completely arbitrary. Moreover, the lack of evidence for a 
very high rotation velocity may mean that the increase in 
lifetime due to rotation may be much smaller than 25 per 
cent in this case, implying small effects from extra mixing 
in the stellar interior. Even if we ignore these considera- 
tions, the corrections just barely make the ages and fiight 
times consistent within (the equivalent to) 2 a. Since this 
star was not corrected for the gravitational darkening effect, 
given the small rotation velocity, it is possible that the dif- 
ference between age and flight time has been overestimated. 
If the correction is applied then the flight time increases to 
46.2tg;g Myr and the age at most to 22.5iio Myr (assuming 
a low metallicity and adding the effect of rotation). Once 
again, we also compared the age determinat i on wit h the one 
obtained form the tracks by iBressan et al.l (|l993h but this 
determination (11.9 Myr) is very similar to the original one. 
We note again that the difference in this situation is still 
significant within the errors. 

HIP 77131 This star was classified as being on the main 
sequence based on the very high (projected) velocity of 
250 kms~^. The a bundan c es of the elements N, C, and Si 
were measured by iMartinI (|2004l ) and, although the abun- 



dances of nitrogen and carbon are consistent with a main 
sequence status, silicon appears to be overabundant. As in 
the case of EC 20252-3137, the assumption of low metallicity 
is not justified in this instance, but even if it were the fiight 
time is still inconsistent with the evolutionary age as can 
be seen in Table O Since this star was not corrected for the 
gravitational darkening effect, as explained in Section 13.31 
once again the difference between age and fiight time may 
have been overestimated. The correction for HIP 77131 was 
done assuming, for the purpose of this analysis, 0.1 dex for 
the gravity and 3000 K for the effective temperature, based 
on a rough extrapolation (assuming a constant function) 
from the values corresp onding to the highe st temperature on 
the theoretical model (jPremat et al.ll2005l ). The fiight time 
obtained after the correction was 17.91^;^ Myr, which is stiU 
inconsistent with the corrected age of 6.3^2 Myr (computed 
assuming low metallicity and adjusted for increased lifetime 
due to high rotation, after correction for the gravitational 
darkening effect), even within twice the error bounds. Again 
in this case we compared the original age esti mate with the 
one obtained using the theoretical tracks by iBressan et al.l 
(1993) but this determination (4.2 Myr) is very similar to 
the original one. 

In summary, all three stars appear to be inconsistent 
with ejection from the disc even if we consider effects which 
would increase their lifetimes. In all cases new observations 
to confirm the atmospheric parameters are highly war- 
ranted, in particular in the case of HIP 77131 since it would 
be very important to have spectroscopic estimates. In that 
respect, the most intriguing case is probably the star SB 357, 
since its atmospheric parameters were obtained from high 
resolution spectra and so should be particularly reliable. 
However, since the composition of the star is unknown, it is 
possible that its metallicty is well below solar, in contrast 
to the other two stars (EC 20252-3137 and HIP 77131). 
For this reason SB s357 is probably the strongest case yet 
for a star formed in the halo. Note that the blue strag- 
gler scenario is not discussed because Mass/Luminosity 
ratios, rotation velocities, and chemical abundances of B 
type blue stragglers in open clusters are indistinguish- 
able from "normal" stars ijSchonberner fc Nataiwotzkil 
1994 lAndrievskv. Schonberner fc Drillind |2000| . 



Schonberner. Andrievskv fc Drillind 200lh . However, 



as was mentioned before in this section, the ejection of 
an intact binary with a high ejection velocity should be 
an extremely rare occurrence. The estimated frac tion of 
binaries among runaway stars is of 1 per cent (|Peretsl 
I2OO9I and references therein), which is a factor of 3 lower 
than the observed fraction (3 in 96 stars). 



5.2 Ejection velocities 

Simulations for bo th ejection mechanisms 

l|Leonard fc DuncanI Il990l and iPortegies ZwartI I2OOOD 
predict a two dimensional mass-energy distribution. In 
particular, this distribution should show an inverse relation 
between mass and ejection velocity. It is known that run- 
away stars will usually be observed near t he apex of th eir 
orbits, where the velocity is the lowest, fsee lMartinll2006l ) - 
explaining why the distribution of measured (heliocentric) 
radial velocities in our sample has an average of only 
41 kms~^. Hence, stars with high ejection velocities will 
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typically be further away than stars with lower ejection 
velocities. This fact introduces a bias in the detection of the 
highest ejection velocity stars, i.e. the stars which constitute 
the tail of the ejection velocity distribution, where the two 
ejection mechanisms are expected to differ. 

The mass-ejection velocity distribution for our sample is 
shown on Fig. |S] with the size of the symbols being propor- 
tional to the (vertical) distance to the disc (larger symbols 
meaning higher distances), and the colour indicating the ap- 
parent visual magnitude. The minimum velocity needed for 
a star to reach a height of 1 kpc above the disc is shown 
as a solid line, explaining why no stars are found below 
said line. This is a result of the selection criteria which 
excluded nearby runaway stars. Nevertheless, the distribu- 
tion appears to show a trend of decreasing ejection velocity 
for higher masses as predicted by theory (the trend is very 
weak and depends on the three more massive stars, includ- 
ing one of the problematic stars mentioned in Section 15.11 
the most massive one). This is the reason why the sam- 
ple is dominated by late type B stars, which is potentially 
problematic because it means that the stars that we need 
to observe are also the fainter. We note that the surveys 
used as a source go as deep as = 15 (considering only 
runaway star candidates), however the faintest runaway ob- 
served has a magnitude of ~ 14.5. We also note that the 
trend of increasing distance to the disc (and corresponding 
decrease in brightness) for higher ejection velocities appears 
to be true only up to about 300 kms^^, with the maxi- 
mum ejection velocity being about 400 - 500 kms~^. The 
fact that we find a mix of bright and relatively faint stars, 
in a magnitude range of 12 < V < 14.5, clustered around 
400 - 500 kms~^ reflects a real drop in ejection velocity 
distribution. In fact, this is the highe s t ejection velocity 
predic ted byjLconard & Duncan| (|l990f ) , iGvarama dze et al] 
l|2009ll . andfPrzvbilla ct al, (20oi|). In Fig.Qwc have plotted 
the cumulative distribution function (CDF) of the distribu- 
tion function and the best Maxwellian flt which peaks at 
156 kms~^. It can be seen that the fit is good for veloc- 
ities up to ~ 300km s~^, however there is a suggestion of 
bimodality which could indicate the existence of two dif- 
ferent populations. If we assume that this group of high 
velocity runaway stars (velocities higher than 350 kms~^, 
corresponding to 11 stars above this threshold and 85 be- 
low) corresponds to a different population and remove them 
from the distribution, then the fit to a Maxwellian distribu- 
tion is much better, peaking at 141 kms~^, as can be seen 
on Fig. [S] Hence, the "slow" group appears to be consis- 
tent with the s tandard ejection scenarios. Moreover, since 
iLeonardl (|l993t ) predicts a Maxwellian distribution peaked 
at 50 - 100 kms~^ for the ejection velocity, the observed 
distribution seems to indicate that we could not be missing 
many objects in the tail of the distribution, if we remember 
that we are missing many low velocity objects because we 
have selected only stars high above the galactic plane. Note 
that the selection of high galactic latitudes may induce an 
overestimate of the number of high velocity stars, since the 
observed cone volume increases with distance. Although this 
effect may explain the strong high velocity tail observed, it 
does not explain the apparent bimodality. Possible selection 
effects and completeness of the sample will be discussed in 
the next section. The significance and possible identity of the 
apparent high velocity population will be discussed later. 
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Figure 6. Ejection velocity - Mass distribution. Tlie apparent 
visual magnitude is given by the grayscale, from black to white, 
black being the brightest and white the faintest. The size of the 
circles is proportional to the height above the galactic plane. The 
line indicates the minimum velocity needed to reach a height of 
1 kpc. 




Figure 7. Ejection velocity distribution. The filled line is the 
empirical cumulative distribution function, the dashed line is the 
best fit Maxwellian distribution (peaking at 156 kms~^) , and the 
dotted line is the the predicted Maxwellian distribution l|Leonardl 
Il993^ (peaking at 100 kms"!). 



5.3 Number estimate and completeness 

An estimate of the proportion of runaways in the local pop- 
ulation of early type stars is important because it provides 
a quantitative measure of the accuracy of the ejection ve- 
locity result, and an important constraint on the theoret- 
ical models. The number of early type stars in the solar 
neighbourhood was estimated from the Hipparcos Catalogue 
(Ivan LeeuwenllioOTl ). by applying a colour cut of B — V < 
(corresponding to a te mperature larger than 9500 K, see 
iNapiwotzki et al.l Il993l ). and a parallax > 5 mas, corre- 
sponding to a distance of 200 pc. The derived surface number 
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Figure 8. Ejection velocity distribution for the sub-sample of 
stars with ejection velocities lower than 350 km s~^ . The filled line 
is the empirical cumulative distribution function, the dashed line 
is the best fit Maxwellian distribution (peaking at 141 kms~^), 
and the dotted line is the the predicted Maxwellian distribution 
(|Leonardlll993 ) (peaking at 100 kms^^). 

density of early type stars in the solar neighbourhood was 
1.1 X lO^kpc"^. The estimate of the local number density 
of runaway stars requires an estimate of the completeness of 
our sample. 

There are a few different effects influencing the com- 
pleteness of the sample: the height above the galactic plane 
and the lifetime on the main sequence are linked to the ejec- 
tion velocity, and to the brightness of the stars. It is im- 
possible to get a firm grasp on the influence of these selec- 
tion effects without doing extensive modeling. Nevertheless 
a rough estimate of completeness is useful for the reasons 
described. The brightness is directly related with the dis- 
tance, which can be decomposed in a radial component and a 
perpendicular component (height above the galactic plane). 
Hence, for a given radius we will start losing stars when we 
go above a certain height. If we assume that the ejection of 
runaway stars is an isotropic process, then the distribution 
of the radial component (projection on the galactic plane) of 
the runaway stars contained inside a cylinder with a given 
radius should be the same as the distribution of bright young 
stars in the disc, inside a circle with the same radius. For 
this reason, we have estimated completeness by computing 
the distribution of the radial component within a circle of 
radius (distance on the galactic plane) r, and comparing it 
with a realistic simulated isotropic distribution of stars in 
the disc, scaled to have the same frequency as the sample at 
radius r. We then varied r until we found a good fit between 
both distributions on the interval - r. Using this method 
we estimated the sample to be complete up to a radius of 
~ 1 kpc. This means that up to a radial component of 1 kpc 
we are complete in the height component of the distance. 
Note that this is a volume limited sample, since we have se- 
lected only stars above ~ 280 pc. This lower limit in height 
coupled with the lower limit of 20° (absolute value) in galac- 
tic latitude explains why we are complete up to a radius of 
1 kpc, which translates into a height of ~ 360 pc for this 
galactic latitude. This is the effect mentioned in Section r5.2l 



which increases the relative number of stars at greater dis- 
tances from the galactic plane or, equivalently the number 
of stars with high ejection velocities. 

From the volume limited complete sample of runaway 
stars we have estimated the local surface number density to 
be 8.6kpc~^. This co rresponds to 27 stars inside a cylinder 
with radius of 1 kpc. iMdzinarishvili fc Chargeishvilil (|2005l ) 
compiled a list of 61 runaway star candidates inside a cylin- 
der of 2 kpc around the Sun using Hipparcos proper motion 
data. Since Hipparcos data is complete only up to V ~ 9 
we know this sample is highly incomplete, and will con- 
sist mostly of nearby stars, however this means that it is 
a good complement to our sample. If we combine the two 
samples we find that the local number density is 13.4 kpc~^. 
Note th at the true value is probably smaller s ince the stars 
from the IMdzinarishvili fc ChargeishvUil l|2005h sample have 
not been confirmed as real runaway stars, nevertheless this 
number is a good estimate of the upper limit. It is inter- 
esti ng to note that ou r estimate is not far from the one 
by [Leonard fc DuncanI ([l990), who estimates a value of 
~ lOkpc-^ 

We have then extrapolated the local densities of run- 
away stars and early type main sequence stars to their total 
number in the galaxy, using the simulated distribution of 
young stars to compute the proportion of the local samples 
relative to the whole galaxy. We estimate the population of 
runaway stars to be ~ 0.1 per cent of the total number of 
young, early type main sequence stars, which appears to be 
a reasonable number, if we assume that the number density 
decreases for later types. This is consistent with estimates 
of space frequencies of runaway stars: 3 per cent fo r O type 
and 4 per cent for the range BO to B5 (|Stondll99ll ). 

5.4 Comparison with other studies 

In this section we divide the literature overlapping with our 
work in four groups, noting that we are interested in com- 
paring results regarding the kinematical analysis: 

(i) the papers dealing with the Palomar-Green sam- 
ple, here abbrev iated to PG papers l|Rolleston et al.lll999l : 
iLvnn et al.ll2004l ): 

(ii) the papers dealing with the Edinburgh-Cape sam- 
ple, here abbreviated to EC papers l|Rolleston et al.l [19971 : 
iMagee et ai]|200ll : iLvnn et al.ll2004l ): 

(iii) the pap ers deal ing with the Hip p arcos sample 
ll Allen fc Kinma n 2004; IConlon et al] Il990l : iMartinI |2004 
l2Q0^; 

(iv) the rem ainin g papers (|Conlon et al.l Il992l : 
iRamspeck et al.ll200ll : lBehj|2003bD . 

The studies focused on the Hipparcos sample are by them- 
selves unsuitable to study the ejection velocity distribution, 
given that we are interested in the tail of the aforemen- 
tioned distribution, where the difference between the BEM 
and DEM should be more evident. The Hipparcos brightness 
limit (see Table |4} is crippling if we want to have an esti- 
mate of the ejection velocity for the more distant, therefore 
potentially faster, objects. This means that although proper 
motions are used in th e kin ematica l analy sis performed by 
lAUen fc KinmanI (|2004l l and lMartinI (|2006h . these studies do 
not give a full picture regarding the nature of the ejection 
velocity distribution. 
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Table 4. Limiting magnitudes and completeness estimates. The 
completeness estimates have a dependency on galactic latitude 
and spectral type. 



Survey 


limiting 


completeness 




magnitudes 


range 


Hipparcos 


V ~ 12.4 


7-3 < V < 9 


Edinburgh-Cape (EC) 


V = 15 


11 < y < 15* 


Palomar-Green (PC) 


n<V < 14.5 


13 ^ BpG ^ 14.6** 



** — Photographic magnitudes BpQ given for PG survey 



On the other hand, the studies dealing with the EC 
and PG surveys suffer from the fact that they explore only 
a narrow range of brightnesses (see Table |4} and proper 
motion information was not used for the kinematical anal- 
ysis of the studied objects. In the EC series of papers, the 
kinematical analysis is incomplete in the sense that proper 
motion data was not available for the stars EC 20011-5005 
and EC 20105-2944. Similarly in the PG series of papers, the 
kinematical analysis lacked proper motion data for five stars 
(40 per cent). The other studies are not systematic, consist- 
ing of observations of only a sm a ll num ber o f objects. The 
stars from the studies by iBehij (|2003lJ ) and IConlon et al.l 
cover a range in magnitude that is barely higher 
than the limiting magnitude from the Hipparcos sample 
and they do not include proper motion data in the kine- 
matical analysis . On th e other hand, the sample studied by 
iRamspeck et al.l (|200ll ) covers the range 10 < 1^ < 14 in 
brightness and does include proper motion data in the kine- 
matical analysis performed for all but four stars: PHL 159, 
PG 1511-^467, SB 357, and HS 1914-H7 139. The star SB 357 
is also analysed bv lConlon et al.l (|l992l ). 

By combining these samples we were able to better con- 
strain the ejection velocity distribution, because we cover 
the near objects (Hipparcos sample) and distant objects 
(the other samples), corresponding to a range in distance 
of 0.57 - 35 kpc. The fact that our kinematical analysis in- 
cludes proper motion data for all stars also constitutes an 
improvement over previous analysis, because it has permit- 
ted us to compute the orbits with all information necessary 
(position and velocity). In relation to the Hipparcos sam- 
ple nothing new has emerged from our study, which is not 
surprising given the fact that the data from Hipparcos was 
already available and it has good quality. On the other hand, 
we have new results on stars from the other samples. All the 
stars mentioned earlier in this section lacking proper motion 
data on previous studies have now more precise estimates of 
the flight times and ejection velocities. The comparison of 
the flight times with the estimated evolutionary ages con- 
firms that most stars are consistent with the runaway hy- 
pothesis, whereas the computed ejection velocities suggest 
a limit of 400 - 500 kms~^ for the ejec tion. This is consis- 
tent w ith t he theoretical pr edictions of ^L eonard &: DuncanI 
l|l990l) and lGvaramadze et al . (2009). The stars SB 357 and 
HS 1914-1-7139 were estimated t o have flight ti mes l arger 
than their evolu tionary ages by IConlon et al.l (|l992l ) and 
|2001), in the case of SB 357, and by 
20011 ) in the case of HS 1914-^7139. The 
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inclusion of proper motion data in the kinematical analysis 
permits us to confirm that SB 357 appears to be incompati- 
ble with ejection from the disc, whereas HS 1914-1-7139 now 
appears to have a flight time compatible with an ejection 
from the disc scenario. It is important to mention that the 
star HS 1914-1-7139 was also studied bv lHeber et al] (|l995l ) 
who confirmed a fiight time much higher than the estimated 
age, based on an analysis lacking prop er motion data. Sim i- 
larly, the star BD -2 3766, analysed bv lConlon et al. I (|l992t ). 
who found a fiight time higher than the estimated age of the 
star, appears to be compatible with ejection from the disc 
after the inclusion of proper motion data in the analysis. 

5.5 The link with hypervelocity stars 

Hypervelocity stars are stars traveling on unbound orbits, 
generally believed to be ejected by the Supermassive Blac k 
Hole in the centre of the Milky Way (iBrown et al.ll2007bl ). 
the so-called Hills mechanism l|Hillsl Il988l ). The apparent 
maximum ejection velocity in our sample is very close to the 
escape velocity in the solar neighbourhood (555.3 kms~^ 
for the adopted gravitational potential). This means that 
the group of stars with the highest ejection velocities (the 
high velocity group discussed on Section 15.21 hereafter ab- 
breviated to high velocity group) are near the regime of hy- 
pervelocity stars. Note that, depending on the direction of 
ejection, some of them could - in principle - exceed the es- 
cape velocity. It is however unclear if they constitute a low 
velocity tail of the hypervelocity population, or if they are 
simply extreme cases of runaway stars. It could also be the 
case that the distinction between hypervelocity and runaway 
stars is artificial and they all should be considered part of 
the same population, the former being just an extreme case 
of the latter. 

There are at least two cases of hypervelocity stars which 
were shown to have be en ejected far from the galactic centre. 
These are HD 271791 llHelaer et al.ll2008l an d lPrzvbiUa et al.l 
[2008) and HIP 60350 (|lrrgang et al.ll201ol ). having ejection 
velocities of ~ 400 kms~^ and 379 kms~^, respectively. 
These velocities correspond to the observed limit in our 
study, as discussed in Section 15.21 which makes the case 
for the existence of a real limit even stronger. Interestingly, 
from the stars in our sample with ejection velocities greater 
than 350 kms~^, only two have points of origin in the disc 
compatible with, although not necessarily coinciding with, 
the centre of the galaxy, HIP 105912 and EC 19596-5356. 
A galactic centre origin can be ruled out for the other nine 
cases. Our results, combined with the results for HD 271791 
and HIP 60350, imply that the high velocity group of run- 
away stars could not, in its entirety, have been ejected by 
the Hills mechanism. 

Both the DEM and the BEM scenarios are able to pro- 
duce stars with the requ ired velocities of ~ 400 km s~^. The 
mechanism proposed bv lGvaramadze et al.l (|2009l ). which is 
a variation of the classical DEM where the dynamical inter- 
action is between a binary and a very massive star (mass 
greater than 50 M0) instead of another binary, can also 
produce bot h runaway and hyp ervelocity stars. In fact it 
is argued by ICvaramadz^ l|2009l ) that the star HD 271791 
achieved its high velocity due to the DEM. However, the 
predicted velocity distribution is continuous in DEM mod- 
els, with no existence of a natural velocity limit (with the 
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excep tion of the 1400 kms~^ limit established by iLeonardI 
ligQll ). These predictions are at odds with the observed dis- 
tribution, where both a limit velocity and a gap in the dis- 
tribution are apparent. A var i ation of the classical BEM as 
described by iPrzvbilla et al.l (|2008l ') has been proposed to 
explain the ejection of both HD 271791 and HIP 60350. In 
this case a very close binary, composed of a massive star 
(mass > 40 Mq) and a B type main sequence, undergoes 
a common envelope phase until the primary loses the enve- 
lope due to the deposition of the spiraling secondary's or- 
bital energy, stopping the merger. The system will then be 
composed of a Wolf-Rayet star and a B type main sequence, 
with a very small separation between the two components, 
which can be as small as 13 — 15 Kq . When the Wolf-Rayet 
ends as a supernova, the B type companion is released with 
a velocity close to its orbital velocity. This scenario estab- 
lishes then a natural limit for the ejection velocity of around 
400 kms~^, close to the observed limit. Moreover, the small 
separation achieved during the common envelope phase, cou- 
pled with the fact that a Wolf-Rayet star is more compact 
than a main sequence star with the same mass, would ex- 
plain the observed gap in the ejection velocity distribution. 
This follows because two different populations would corre- 
spond to different types of binaries: a binary with a very 
massive Wolf-Rayet star, evolving into a common envelope 
phase, as previously de scribed, or the more traditional sce- 
nario (e.g. IPortegies Z wart , 2000. ') where the binary evolves 
either as a detached or semi-detached binary. Note that the 
shorter separation in the latter case implies higher ejection 
velocities. 

On the other hand, it could be argued that the tradi- 
tional mechanisms can only accelerate stars up to velocities 
of about 200 kms~^ as has been tra diti onally believed (see 
for e xample iHoogerwerf et all I2OOII and iGvaramadze et al] 
I2OO9I ). what would mean that the stars with high ejection 
velocities in our sample were ejected by a different mech- 
anism. This mechanism would be a third mechanism re- 
sponsible for the ejection of hypervelocity stars far from 
the galactic centre, with the stars with high ejection ve- 
locities corresponding to the low velocity tail of the ejec- 
tion velocity distribution produced by this mechanism. A 
ca ndidate would be for example the m echanism proposed 
bv iGualandris fc Portegies Zwar"3 l|2007l ). which consists in 
a dynamical encounter between a binary containing the 
ejected star and an intermediate-mass black hole in the cen- 
tre of dense cluster. The fact that a gap in the ejection ve- 
locity distribution is observed suggests that they may cor- 
respond to two different populations, however only a better 
understanding of all selection effects would permit to draw 
a firm conclusion. This hypothesis suffers from the fact that 
direct evidence of the existence of intermediate-mass black 
holes is still lacking. 

In summary, everything considered, the evidence sug- 
gests that there is no reason to separate runaway and hyper- 
velocity stars in two separate populations. All mechanisms 
have to be considered together to understand the whole pic- 
ture, including the BEM, the DEM, the Hills mechanism, 
and possibly the intermediate-mass black hole interaction 
mechanism. The most extreme ejection velocities are better 
explained by the Hills mechanism when they occur in the 
galactic centre, and by the BEM when they occur far from 
the galactic centre. 



6 CONCLUSIONS 

With this work we wanted to achieve a better understanding 
of the mechanisms behind the runaway stars phenomenon. 
We aimed at producing an estimate of the ejection velocity 
distribution to serve as a constraint on the theoretical mod- 
els. Stars far away from the galactic plane are ideal for this 
purpose because they have been ejected with higher veloci- 
ties than most runaway stars still traveling in the disc. This 
is important since the models differ mostly in the high veloc- 
ity of their ejection velocity distributions. We have greatly 
improved on previous analyses by: combining all the pre- 
viously existing samples into a homogeneized one, includ- 
ing hitherto unavailable proper motion measurements in the 
kinematical analysis, doing a systematic treatment of the er- 
ror propagation, and correcting for the effect of gravitational 
darkening. The main conclusions are: 

(i) a small sub-sample of three stars (SB 357, EC 20252- 
3137 and HIP 77131) cannot be directly explained by any 
of the ejection mechanisms, as they have flight times much 
larger than their evolutionary lifetimes. This is made clear 
by the ex ample of the sta r SB 3 57 that was already sus- 
pected bv iRamspeck et al.l l|200ll '). but whose case is made 
here more convincingly because we have included the proper 
motion of the star in the kinematical analysis. An indepen- 
dent confirmation of the atmospheric parameters is highly 
warranted. In contrast, some stars which appeared to be 
incompatible with ejection from the disc, are shown to be 
compatible once proper motions arc included in the analy - 
sis. Note, for example, the cases of PHL 346 llConlonlll99j ) 
and HS 1914-^7139 (|Ramspeck et al.ll200ll ). 

(ii) we determined the ejection velocity distribution, 
which shows evidence for bimodality and a limiting ve- 
locity of 400 - 500 km s~^ The BEM, as presented by 
IPrzvbilia et all l|2008h . seems to be enough to produce the 
observed ejection velocity distribution. In this scenario, 
the limiting velocity seems to correspond to a real limit 
which arises naturally. The bimodality in the distribution 
could also be a natural consequence of the BEM. How- 
ever, a contribut i on fro m the DEM was demonstrated by 
IHoogerwerf et al] (|200ll ). Note that the observed group of 
high velocity stars makes the observed distribution devi- 
ate from the expected Maxwellian. More extensive modeling 
would be needed to assess if the mechanisms are able to pro- 
duce the observed numbers. 

(iii) the stars with the highest ejection velocities may be 
instead hypervelocity stars, at least in the general sense of 
corresponding to the lower velocity tail of a population of 
hypervelocity stars. Alternatively, hypervelocity stars can 
be seen as more extreme cases of runaway stars. There is 
at least some evidence for two populations given by the bi- 
modality in the distribution. We have also determined that 
most of these high velocity stars are inconsistent with an ori- 
gin close to the galactic centre. These results are consistent 
with the two hypervelocity stars which had their ejection 
velocities and places o f origin in the gal acti c plane deter- 
mine d in other studies (|Heber et al.ll2008 and'lrrgan g et akl 
i2010[). Since the BEM appears to be able to explain the 
features of the ejection velocity distribution, we conclude 
that it is probably the mechanism producing hypervelocity 
stars. It would be interesting to do a kinematical analysis of 
the remaining sample of hypervelocity stars and see if the 
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ejection velocities and places of origin are still consistent 
with this scenario, however this is dependent on the delicate 
measurement of proper motions of very distant stars. The 
launch of the Gaia mission in 2012 will hopefully permit us 
to approach this problem. 
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APPENDIX A: THE GALACTIC POTENTIAL 

The d ynamical mass model created by lAllen & SantillanI 
(jl993), hereafter called the default model, has been used 
as a basis for several dynamical studies. It is constructed 
as a sum of three independent components corresponding 
roughly to the bulge, disc, and halo populations, with the 
dark mass contribution included in the Halo component. The 
model successfully reproduces the galaxy rotation curve and 
its main advantage is the fact that the functional form of 
the gravitational potential for the three components is quite 
simple, being also fully analytical, continuous, and with con- 
tinuous derivatives everywhere. However the fact that it was 
created in 1991, in the pre-Hipparcos age, means it has some 
problems which we have addressed by updating it so it agrees 
with more recent results on the observable quantities. The 
problems we identified were: the disc component scale-length 
was too high, the local density was too high, as was the in- 
tegrated surface density. 

The disc scalelength used on the default model has the 
the value of 02 — 5.3178. this value is a factor of 2 higher 
than recent estimat es obtained from analysis of tracers of 
young populations. iDrimmel & Spergell (|200ll ) obtained a 
value of a2 = 2.25 from infrare d observat i ons tr acing the 
galactic dust in the disc, whereas ISale et al] l|2010l ) obtained 
a value of 02 = 3.0 by counting A-type stars in the outer 
disc. Changing the scale-length to the appropriate value in 
the default model produces ans excess of mass in the centre 
of the galaxy. 

As a means of reducing the mass concentrated near the 
centre of the galaxy we introduced a second disc component, 
wi th negative mass, an alogous to the "hole" component used 
bv lRobin etHI (|2003l ') for their density model. In this case 
this component can also be thought of as a "hole" in the 
original disc, at least in the sense that the disc density is 
greatly reduced in the central regions. This new component 
has a functional form equal to the "real" disc component, 
hereafter first disc component, only with a negative sign, 
which ensures that we keep the simplicity and good math- 
ematical properties of the default model. We designate this 
new model with four components by "updated model" , and 
the new disc component by "second disc component". 

To distinguish the parameters of the first and second 
disc components a superscript with the number 1 or 2, re- 
spectively, is used. The parameters used to constrain the 
model are summarized in Table \KT\ The procedure adopted 
to construct the updated model consisted of the following 
steps: 
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Table Al. Physical quantities that were used to constraint the 
parameters of the models, compare d with the measur e d val - 
ues. The local density value is from iHolmbere fc FlvnrJ 1I2OOOI ) 
and the circular ve l ocity at the position of the Sun is from 
iBinnev &: Tremainj 1I2OO8I ). Values for two different distances rg 
are given to facilitate the comparison between the models,- the 
distance to the Sun is assumed to be 8 kpc in the updated model, 
but 8.5 kpc in the default model. 



Table A3. Comparison of physical quantities computed from the 
model with measu red values. The Port consta nts and the escape 
velocity are from Binncv & Trcmainc ('200^), the local surface 
densi ty is from .Holmbcrg fc Flymi (2004), and the total mass 
from lXue et al. Values for two different distances ro are 

given to facilitate the comparison between the models - the dis- 
tance to the Sun is assumed to be 8 kpc in the updated model, 
but 8.5 kpc in the default model. 





default 




updated 


reference 




model 




model 


values 


ro (kpc) 


8.5 


8.0 


8.5 8.0 


8.0 


po (M0PC-3) 


0.15 


0.16 


0.096 0.112 


0.102 ± 0.010 


v'^ (kms^-^) 


220 


220 


218 217 


220 



Table A2. Updated model parameters and total masses derived 
for each component. 



default 



updated 



reference 





model 




model 




values 


ro (kpc) 


8.5 


8.0 


8.5 


8.0 


8.0 


So (Mq pc-2) 


83 


93 


62 


71 


74 ±6 


A (kms~l kpc~^) 


13.0 


13.5 


13.6 


14.5 


14.8 ±0.8 


B (kms-l kpc^i) 


-12.9 


-13.9 


-11.9 


-12.7 


-12.4 ±0.6 


v° (kms-i) 


536 


541 


550 


555 


550 ± 50 


Mtot (Mq) 


9.00 X 10" 


1.00 X 10^2 


~ 1.00 X 10^2 



Component Mass (Mq) Parameters 

Central mass 1.41 X 10^" Mi = 606.0 

bi = 0.3873 



First disc 



Second disc 



7.14 X 10^ 



= 3079.0 
al = 3.0 
fei = 0.25 



M| = -307.9 
-7.14 X 10^ a'j = 1.5 
= 0.25 



Halo 



9.00 X 10^ 



Ms = 4761.1 
as = 12.0 



• the first disc component scalelength, was set to 3 kpc, 
for tlie reason explained before. 

• tiie second disc component was established first, by 
choosing suitable parameters to satisfy two conditions: re- 
duce the disc mass accumulated in the centre by as much as 
possible, while keeping the rotation curve as flat as possible 
in the region between 5-20 kpc. The vertical scale height 
was kept fixed and equal to the one in the first disc compo- 
nent. The parameter Ml was chosen to be 10 per cent the 
parameter M2. 

• keeping the second disc and bulge components fixed, a 
first attempt at fitting a flat rotation curve over the interval 
5-20 kpc, and slowly decreasing until 100 kpc, by varying 
the first disc mass (parameter M2), and the halo compo- 
nent mass and power law exponent (parameters Ms and 7, 
respectively). The estimate obtained for the parameter Ms 
resulting from this fit was taken as definitive. 

• finally, estimates of the 7 and AI2 (and so also of M|) 
parameters were obtained by simultaneously fitting a flat 
curve to the interval 5-20 kpc, imposing the observed local 
circular velocity, and the observed local density. 

A summary of the parameters obtained through the de- 
scribed fitting procedure for the updated model can be seen 
in Table IA21 together with the mass of each component. 
The quality of the model can be judged from Fig. lAll and 




Figure Al. Plot of the circular velocity as a function of ra- 
dius over the interval 0-20 kpc. The dashed line co rresponds to 
the de fault model, the dotted line correspond to the lRobin et al.l 
l|2003l ) model, the solid line corres ponds to the updated m odel, 
the error bars are data points from lAUen fc SantillarJ l ll99ll) . 



from Tableau In Fig.ETIwe can see a comparison b e tween 
the updated model and the default and lRobin et al.l (|2003l ) 
models. The updated model behaves similarly to the others 
in the critical region between 5-20 kpc. Inside 5 kpc the 
behaviour is similar to the default model, which should not 
be surprising, since it is based on it. Note in particular that 
the agreement with the observed data points is quite good, 
although no effort was made to perform an explicit fit. 

Table IXSl verifies that the agreement of quantities com- 
puted from the model - local surface density, Oort's A and 
B constants, escape velocity, and total mass of the galaxy - 
is good, in particular when compared with the performance 
of the default model. Note that the default model was fine- 
tuned assuming a distance of the Sun from the galactic cen- 
tre of 8.5 kpc, whereas we assumed a distance of 8.0 kpc. 
We conclude that the updated model is a better representa- 
tion of the galactic gravitational potential than the default 
model. 
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APPENDIX B: DATA FOR MAIN SEQUENCE 
STARS 

In this Appendix we present the relevant data for the stars 
classified as main sequence. In Tables |BT] and |B2] the val- 
ues used as input to compute the orbits are shown for 
the groups A and B respectively. The proper motions were 
obtained from a combination of the values from UCAC 2 
JZacharias et al.1 120041), USNO-B (iMonet et al.l l2003h. 



Tych o-2 ( H0g et al.l I2OOOI ) and Hipparcos ( van LeeuwerJ 



l200'<t ). SuperCOSMOS, S PM (|Girard et all |2004| ). and 



NPM2 iHanson et all 120031 '). by a weighted_ average follow- 
ing the prescription in Section 2 of iPauli et al.l (|2006l V The 
gravities, temperatures and projected rotation velocities 
were taken from the literature or computed from Stromgren 
uvbyP photometry as indicated. Note that these are the mea- 
sured gravities and temperatures not corrected for the grav- 
itational darkening effect. The radial velocities are heliocen- 
tric correcte d. The masses were d erived from the evolution- 
ary tracks of lSchaller et al.l (|l992l ). See Section [3] for details 
of the procedure. 

In Tables lB3l and lB4l data concerning the present con- 
dition of the stars (present coordinates {X,Y,Z), distances 
and evolutionary ages) and the condition at the moment 
of ejection from the disk (point of ejection in the plane 
{Xq,Yo), flight time and ejection velocities) are shown. All 
coordinates are given in a galactocentric Cartesian right- 
handed frame of reference where the X axis passes through 
the position of the Sun, pointing in the opposite direction, 
and Z points to the north galactic pole. Note that the dis- 
tance from the galactic centre to the Sun was assumed to be 
8 kpc. Full three-dimensional velocities were computed from 
the knowledge of the distances, proper motions and radial 
velocities. The ejection velocities correspond to the orbital 
velocity for Z — 0, corrected for the disk rotation velocity 
component. The evolutionary ages were obtained by in ter- 
polating the evolutionary tracks of lSchaller et al.l (|l992l ). 

The asymmetric error bounds correspond to confidence 
intervals enclosing the equivalent to 1 — cr errors of the er- 
ror distribution computed using the Monte Carlo procedure 
detailed in Section [l] 
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Table BI. Input data for stars belonging to group A. The data shown are: the radial velocity, the proper motions in right ascension 
and declination directions, the surface gravity, temperature, mass and projected rotation velocity. 



Name 


r.v. 


p.m. a cos S 


p.m. (5 


logs 




T'oft 




Mass 


V sin i 




(km s ^ ) 


(mas yr ^ ) 


(mas yr ^) 


(cm s ) 




(K) 






(km s 1 ) 


EC 19337-6743 


-17 ± 10 (14) 


6.24 ± 0.50 


-5.64 ± 0.61 


3.90 ± 0.20 


(14) 


10000 ± 2000 


(14) 


2.8i:°? 


200 (14) 


EC 05515-6231 


-12 ± 10 (14) 


-1.36 ± 0.86 


6.49 ± 0.85 


4.20 ± 0.20 


(14) 


14000 ± 2000 


(14) 




65 (14) 


EC 20089-5659 


-24 ± 10 (14) 


3.61 ± 0.70 


-7.97 ± 0.72 


3.90 ± 0.20 


(14) 


14000 ± 2000 


(14) 


4 4+11 


90 (14) 


EC 03462-5813 


39 ± 10 (19) 


3.33 ± 0.64 


21.77 ± 0.65 


4.20 ± 0.20 


(19) 


13500 ± 675 ( 


L9) 




200 (19) 


EC 05582-5816 


81 ± 10 (19) 


9.85 ± 0.91 


10.10 ± 0.84 


4.00 ± 0.20 


(19) 


17000 ± 850 ( 


L9) 


5 6 + °-^ 


170 (19) 


EC 10087-1411 


105 ± 10 (19) 


-9.34 ± 0.76 


-5.92 ± 0.80 


4.10 ± 0.20 


(19) 


14500 ± 725 ( 


L9) 


4 2 + 0-5 
— 5 


180 (19) 


EC 05490-4510 


32 ± 10 (19) 


-1.75 ± 1.34 


2.83 ± 1.37 


4.20 ± 0.20 


(19) 


17000 ± 850 ( 


19) 


5 1 + "-^ 
^ — 5 


30 (19) 


EC 10549-2953 


-17 ± 10 (19) 


-8.81 ± 1.09 


1.03 ± 1.07 


4.10 ± 0.20 


(19) 


14000 ± 700 ( 


19) 


4 0+0-5 


200 (19) 


EC 19071-7643 


-20 ± 10 (14) 


16.66 ± 0.65 


-20.15 ± 0.64 


4.20 ± 0.20 


(14) 


21500 ± 2000 


(14) 




30 (14) 


EC 19476-4109 


-19 ± 10 (14) 


-3.70 ± 0.61 


-4.51 ± 0.62 


4.00 ± 0.20 


(14) 


17000 ± 2000 


(14) 




120 (14) 


EC 20292-2414 


15 ± 10 (19) 


1.36 ± 0.78 


-11.32 ± 0.85 


4.10 ± 0.20 


(19) 


25000 ± 1250 


(19) 


10.3+1-^ 


160 (19) 


EC 13139-1851 


23 ± 10 (19) 


-9.15 ± 1.39 


-3.80 ± 1.52 


4.20 ± 0.20 


(19) 


17000 ± 850 ( 


19) 


+ 0.6 
— 6 


50 (19) 


EC 20140-6935 


17 ± 10 (19) 


8.81 ± 0.51 


1.79 ± 0.62 


3.70 ± 0.20 


(19) 


20500 ± 1025 


(19) 




65 (19) 


EC 09414-1325 


71 ± 10 (19) 


-2.36 ± 1.03 


-1.47 ± 1.10 


4.10 ± 0.20 


(19) 


14000 ± 700 ( 


19) 


3.9l°o1 


260 (19) 


EC 03240-6229 


-12 ± 10 (19) 


3.27 ± 1.23 


4.28 ± 1.28 


3.80 ± 0.20 


(19) 


12000 ± 600 ( 


19) 




165 (19) 


EC 20153-6731 


-40 ± 10 (19) 


5.07 ± 1.07 


-8.76 ± 1.00 


3.80 ± 0.20 


(19) 


14500 ± 725 ( 


19) 


4 8+°-^ 

— 6 


120 (19) 


EC 19586-3823 


-102 ± 10 (14) 


-0.13 ± 1.04 


2.56 ± 1.03 


3.90 ± 0.20 


(14) 


18500 ± 2000 


(14) 




150 (14) 


EC 06387-8045 


60 ± 10 (19) 


-0.99 ± 0.92 


6.62 ± 0.94 


3.90 ± 0.20 


(19) 


22500 ± 1125 


(19) 


9.5+_\i 


190 (19) 


PC 1533-1-467 


33 ± 6 (18) 


-11.45 ± 0.53 


11.41 ± 0.54 


4.09 ± 0.10 


(18) 


18500 ± 925 ( 


18) 


6 0+°-^ 


215 (18) 


EC 20252-3137 


26 ± 10 (14) 


-6.08 ± 0.89 


-3.59 ± 0.91 


4.00 ± 0.20 


(14) 


23000 ± 2000 


(14) 


9 2+^-^ 
— 1 5 


60 (14) 


EC 05438-4741 


53 ± 10 (19) 


-0.35 ± 1.07 


-0.94 ± 1.53 


4.10 ± 0.20 


(19) 


13500 ± 675 ( 


19) 


„ O + 0.5 
3.0_Q 4 


30 (19) 


EC 20104-2944 


135 ± 10 (14) 


-0.45 ± 1.70 


-1.40 ± 1.98 


4.20 ± 0.20 


(14) 


15000 ± 2000 


(14) 


4.2lJ-°9 


50 (14) 


EC 05515-6107 


89 ± 10 (19) 


-5.82 ± 1.65 


11.57 ± 1.43 


4.00 ± 0.20 


(19) 


22000 ± 1100 


(19) 




290 (19) 


PC 1205-1-228 


153 ± 4 (20) 


-15.44 ± 0.55 


0.05 ± 0.48 


4.10 ± 0.20 


(20) 


16600 ± 1000 


(20) 


5 2+°-^ 

— 6 


165 (20) 


EC 23169-2235 


82 ± 10 (19) 


4.24 ± 1.25 


4.35 ± 1.18 


4.40 ± 0.20 


(19) 


15000 ± 750 ( 


19) 


3 9+0-^ 

— 4 


140 (19) 


EC 09452-1403 


236 ± 10 (19) 


-3.55 ± 1.54 


-0.37 ± 1.71 


4.30 ± 0.20 


(19) 


14000 ± 700 ( 


19) 




70 (19) 


PC 2345-1-241 


80 ± 3 (20) 


-2.22 ± 0.74 


-3.09 ± 0.67 


4.20 ± 0.20 


(20) 


18800 ± 1000 


(20) 




54 (20) 


PHL 159 


88 ± 3 (18) 


-3.02 ± 0.87 


-9.38 ± 0.69 


3.59 ± 0.10 


(18) 


18500 ± 925 ( 


L8) 


7 q+OS 
— 8 


21 (18) 


EC 10500-1358 


99 ± 10 (19) 


-4.00 ± 0.93 


-2.17 ± 1.01 


3.80 ± 0.20 


(19) 


15000 ± 750 ( 


19) 


5 1+°'^ 


100 (19) 


PC 1511-1-367 


102 ± 11 (18) 


-3.64 ± 0.88 


-13.10 ± 0.87 


4.15 ± 0.10 


(18) 


16100 ± 805 ( 


18) 


4 7+0'-i 
— 4 


77 (18) 


EC 20011-5005 


-171 ± 10 (14) 


2.06 ± 1.61 


-3.96 ± 1.83 


4.20 ± 0.20 


(14) 


17000 ± 2000 


(14) 




30 (14) 


BD -15 115 


93 ± 4 (18) 


7.04 ± 1.00 


0.46 ± 0.76 


3.81 ± 0.10 


(18) 


20100 ± 1005 


(18) 




35 (18) 


PC 2219-1-094 


-24 ± 9 (18) 


-3.53 ± 1.17 


-8.43 ± 1.16 


3.58 ± 0.10 


(18) 


19500 ± 975 ( 


18) 


o ,+0.8 


225 (20) 


HS 1914+7139 


-39 ± 10 (6) 


-0.66 ± 1.86 


-0.14 ± 2.60 


3.90 ± 0.10 


(18) 


17600 ± 880 ( 


18) 


6 1+°-^ 
°'^-0.5 


250 (18) 


PC 0855-1-294 


65 ± 4 (20) 


5.86 ± 0.79 


-3.98 ± 0.82 


3.80 ± 0.20 


(20) 


18100 ± 1000 


(20) 


6-8Io.8 


110 (20) 


PC 0955-1-291 


76 ± 14 (20) 


-7.84 ± 0.96 


-4.27 ± 1.14 


3.80 ± 0.20 


(20) 


13600 ± 1000 


(20) 


4 3+°0 


190 (20) 


PC 2229-H099 


-22 ± 5 (20) 


-1.22 ± 1.85 


-0.98 ± 2.35 


4.00 ± 0.20 


(20) 


17600 ± 1000 


(20) 


5 9+° ** 


16 (20) 


PC 1610-1-239 


91 ± 10 (18) 


-0.20 ± 2.80 


-5.20 ± 1.80 


3.72 ± 0.10 


(18) 


15500 ± 775 ( 


18) 




75 (18) 


PC 2111-1-023 


-153 ± 11 (12) 


-6.98 ± 1.76 


-4.07 ± 0.97 


3.40 ± 0.10 


(12) 


14000 ± 700 ( 


12) 


5 2+°-^ 


140 (12) 


PG 0009-1-036 


142 ± 18 (20) 


-0.99 ± 1.88 


1.43 ± 2.34 


3.60 ± 0.20 


(20) 


14800 ± 1000 


(20) 




350 (20) 


PC 0122-f214 


26 ± 5 (18) 


-1.10 ± 0.58 


-3.23 ± 0.84 


3.86 ± 0.10 


(18) 


18300 ± 915 ( 


L8) 




117 (18) 


SB 357 


58 ± 10 (18) 


0.25 ± 1.27 


-2.81 ± 1.57 


3.90 ± 0.10 


(18) 


19700 ± 993 ( 


18) 




180 (18) 


EC 04420 -1908 


205 ± 10 (14) 


-1.27 ± 1.81 


-1.65 ± 2.23 


3.40 ± 0.20 


(14) 


14000 ± 2000 


(14) 




180 (14) 


PG 1332-H137 


148 ± 10 (17) 


-5.99 ± 0.77 


-7.47 ± 0.72 


3.50 ± 0.20 


(12) 


15000 ± 750 ( 


L2) 




140 (12) 


EC 19596-5356 


200 ± 15 (13) 


1.22 ± 1.95 


-0.52 ± 2.47 


3.75 ± 0.20 


(13) 


16500 ± 1500 


(13) 


5 9+1" 


250 (13) 


PHL 346 


63 ± 4 (18) 


4.12 ± 1.23 


-7.94 ± 1.17 


3.58 ± 0.10 


(18) 


20700 ± 1035 


(18) 


9 6+1-^ 


45 (18) 


PG 0914-1-001 


90 ± 20 (20) 


1.15 ± 1.80 


1.64 ± 2.72 


3.10 ± 0.20 


(20) 


12300 ± 1000 


(20) 


c O + 0.8 
°-'^-0.7 


325 (20) 


PG 1209-1-263 


52 ± 11 (12) 


0.40 ± 1.92 


-0.74 ± 2.38 


3.00 ± 0.20 


(12) 


12000 ± 600 ( 


L2) 




70 (12) 



References: ( 1) thi s study, frorn Strom gre n m)fe);/3 photometry; (2) iBehd ll2003"al): f 3) [Conlon et al.l lll989l'): (4 ') 'Conlon et ahl lll992l'): f5' 



lEvand lll967t ): CellHeber et al.l lll995l'l : (71 iKeenan et alJ lll982l): (81 iKeenan fc DuftorJ (Il983l l: 

lll992ll : (111 'Lcvcnha ren fc Leist er' ('200eO; (12) 'Lynn et al.' ('2004'); (13) 'Lynn et al. ('2004'); _ „ ^ 

l|2003ll : (161 iMartin (,20(3i l: (171 [Martin. 1,2006,); (181 .Ramspeck et al.. (^2001, 1: ( 191 .RoUeston et al., (|1997l l: (201 lRolleston"et al., (|1999i l 



(9) [Keenaii^ al. 

(14) M agee et alJ 



3); (lO)lKilkenir 
y); (15) 'Martin' 



I 
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Table B2. Input data for stars belonging to group B. The data shown are: the radial velocity, the proper motions in right ascension and 
declination directions, the surface gravity, temperature, mass and projected rotation velocity. 



Name 


r.v. 


p.m. a cos S 


p.m. (5 


log 9 
















Mass 


V sin i 




(kms^l) 


(mas yr ^ ) 


(masyr ^) 


(cm s 


2 






(K) 






(Mq) 


(km s 1) 


HIP 81153 


146 ± 7 (17) 


-8.64 ± 0.65 


-13.35 ± 0.65 


4 


05 


± 





45 


(16) 


29891 


± 


684 (16) 




80 (16) 


HIP 1241 


± 5 (17) 


5.98 ± 0.63 


-0.43 ± 0.49 


3 


98 


i 





15 


(1) 


11489 


± 


574 (1) 


q 1 +0.3 

— 3 


130 (16) 


HIP 28132 


54 ± 3 (17) 


-3.90 ± 0.65 


8.35 ± 0.71 


4 


00 


i 





15 


(16) 


16948 


± 


140 (16) 


5 5+"'*' 


100 (16) 


HIP 2702 


-26 ± 10 (14) 


-5.23 ± 0.56 


-5.77 ± 0.57 


4 


20 


i 





20 


(14) 


15500 


± 


2000 (14) 


4-4lii, 


30 (14) 


HIP 55051 


20 ± 10 (8) 


28.82 ± 0.50 


20.03 ± 0.41 


4 


06 


± 





25 


(1) 


25518 


± 


1276 (1) 




150 (11) 


HIP 114690 


83 ± 7 (7) 


4.98 ± 0.71 


-12.56 ± 0.51 


3 


50 


+ 





25 


(9) 


20008 


± 


1000 (1) 




238 (7) 


HIP 11844 


-45 ± 6 (17) 


6.13 ± 0.90 


-2.76 ± 0.95 


4 


14 


± 





08 


(16) 


11202 


± 


68 (16) 


2 9+"l 

— 1 


130 (16) 


HIP 16130 


28 ± 6 (17) 


8.65 ± 0.63 


2.59 ± 0.70 


4 


11 


i 





15 


(1) 


14061 


± 


703 (1) 


^^ — 4 


60 (16) 


HIP 58046 


20 ± 2 (17) 


8.12 ± 0.83 


6.73 ± 0.82 


4 


25 


i 





25 


(16) 


13357 


± 


668 (16) 


q c+0.5 
■-^■^ — 4 


230 (16) 


BD +36 2242 


2 ± 4 (2) 


-5.92 ± 0.52 


1.68 ± 0.50 


3 


87 


i 





15 


(1) 


11806 


± 


590 (1) 


^ — 3 


77 (2) 


HIP 98136 


29 ± 6 (17) 


-3.78 ± 0.72 


-10.84 ± 0.56 


3 


42 


± 





20 


(1) 


21570 


± 


1079 (1) 


11.6+11 


140 (16) 


HIP 111396 


-33 ± 6 (17) 


-0.76 ± 0.62 


2.89 ± 0.62 


3 


62 


i 





15 


(1) 


13662 


± 


683 (1) 




35 (16) 


HIP 77131 


-9 ± 5 (5) 


-4.59 ± 0.75 


-1.76 ± 0.71 


3 


86 


i 





25 


(16) 


32089 


± 


536 (16) 


19.5^=5 


250 (16) 


HIP 113735 


-31 ± 10 (9) 


10.12 ± 0.49 


18.42 ± 0.49 


4 


12 


± 





20 


(1) 


21645 


± 


1082 (1) 


7 9+1" 


? 


HIP 59955 


28 ± 10 (17) 


-10.82 ± 0.72 


-2.41 ± 0.59 


3 


66 


i 





15 


(1) 


13860 


± 


693 (1) 


4 7+0-5 
— 4 


225 (16) 


HIP 76161 


-51 ± 7 (15) 


0.61 ± 0.76 


-0.39 ± 0.71 


3 


43 


i 





20 


(1) 


22331 


± 


1117 (1) 


T9 7+3.0 

— 2 


? 


HIP 96130 


-51 ± 6 (17) 


1.32 ± 0.72 


-0.04 ± 0.58 


3 


63 


i 





21 


(16) 


23141 


± 


239 (16) 


ii.stli 


60 (16) 


HIP 61800 


-48 ± 13 (17) 


0.07 ± 0.45 


-5.11 ± 0.43 


3 


24 


i 





15 


(1) 


12247 


± 


612 (1) 




150 (16) 


HIP 114569 


94 ± 5 (17) 


46.35 ± 0.79 


32.73 ± 0.74 


4 


12 


i 





15 


(1) 


18518 


± 


926 (1) 


6 1+°-'* 

— 6 


70 (16) 


HIP 11809 


± 10 (17) 


-18.90 ± 0.88 


-13.46 ± 0.87 


4 


25 


± 





25 


(16) 


13265 


± 


663 (16) 


c, c + O.S 

— 4 


240 (16) 


HIP 16758 


98 ± 10 (14) 


7.37 ± 0.48 


10.04 ± 0.53 


3 


60 







20 


(14) 


23453 


± 


1173 (1) 


12 1 + ^-5 
— 1 8 


90 (14) 


HIP 108215 


-79 ± 7 (7) 


8.23 ± 0.59 


-8.49 ± 0.38 


3 


06 







15 


(1) 


13996 


± 


700 (1) 




169 (7) 


HIP 115347 


23 ± 7 (7) 


-6.30 ± 0.59 


-7.01 ± 0.46 


3 


49 







20 


(1) 


21106 


± 


1055 (1) 


lo.etli 


27 (7) 


HIP 115729 


26 ± 3 (17) 


3.05 ± 0.72 


-1.21 ± 0.61 


3 


46 







15 


(1) 


17966 


± 


898 (1) 




25 (16) 


HIP 3812 


19 ± 10 (14) 


0.80 ± 0.70 


-1.73 ± 0.73 


4 


00 


± 





20 


(14) 


19000 


± 


2000 (14) 




219 (7) 


PG 1530+212 


-7 ± 25 (2) 


4.04 ± 0.52 


6.60 ± 0.53 


4 


00 







25 


(2) 


15000 


± 


500 (2) 


4 6+°-^ 


104 (12) 


HIP 79649 


19 ± 6 (17) 


-0.48 ± 0.57 


-2.97 ± 0.62 


3 


69 







20 


(1) 


21721 


± 


1086 (1) 




90 (16) 


HIP 12320 


24 ± 12 (17) 


4.33 ± 0.84 


-1.25 ± 0.92 


4 


00 







25 


(16) 


13382 


± 


669 (16) 


3 9+" '^ 

— 5 


225 (16) 


HIP 109051 


72 ± 8 (17) 


1.43 ± 0.72 


-0.88 ± 0.74 


3 


85 







20 


(1) 


21021 


± 


1051 (1) 


o o + l.'3 
8.2_j^ 


105 (16) 


HIP 111563 


42 ± 7 (7) 


5.13 ± 0.47 


-17.86 ± 0.36 


3 


11 







20 


(1) 


22545 


± 


1127 (1) 


17.3+^-;! 


105 (7) 


HIP 13800 


-6 ± 2 (2) 


4.69 ± 0.68 


-3.63 ± 0.80 


3 


67 







15 


(1) 


16632 


± 


832 (1) 




12 (2) 


HIP 48394 


14 ± 12 (17) 


3.34 ± 0.93 


-5.69 ± 0.72 


3 


75 


± 





25 


(16) 


20021 


± 


1001 (16) 




225 (16) 


HIP 107027 


117 ± 8 (17) 


-6.18 ± 0.77 


-11.99 ± 0.64 


3 


65 


± 





21 


(16) 


21959 


± 


231 (16) 


10.5+1-^ 


150 (16) 


HIP 1904 


-37 ± 10 (9) 


2.95 ± 0.47 


-6.27 ± 0.49 


3 


03 


± 





15 


(1) 


12894 


± 


645 (1) 




? 


HIP 55461 


72 ± 6 (17) 


-0.97 ± 0.74 


-5.36 ± 0.69 


4 


00 


± 





20 


(3) 


15600 


± 


1000 (3) 


. .., + 7 
^■^-0.6 


140 (16) 


HIP 112790 


15 ± 4 (17) 


-18.34 ± 0.67 


-13.53 ± 0.69 


3 


53 


± 





15 


(1) 


15628 


± 


781 (1) 


5 9+"-" 


70 (16) 


HIP 59067 


35 ± 5 (17) 


-8.47 ± 0.83 


-5.87 ± 0.70 


3 


89 


± 





15 


(16) 


14974 


± 


98 (16) 


4 8+"-" 


70 (16) 


BD +20 3004 


22 ± 14 (2) 


-10.91 ± 0.51 


5.28 ± 0.55 


3 


40 


± 





15 


(1) 


14070 


± 


704 (1) 


5 6+"-** 


105 (2) 


HIP 105912 


-17 ± 7 (7) 


10.16 ± 0.57 


-11.61 ± 0.47 


3 


29 


± 





20 


(1) 


22201 


± 


1110 (1) 


13.6+2-^ 


102 (7) 


HIP 45904 


36 ± 14 (17) 


3.27 ± 0.79 


-9.63 ± 0.60 


3 


48 


± 





20 


(1) 


21016 


± 


1051 (1) 


10.7t.li 


277 (12) 


HIP 37903 


84 ± 10 (17) 


2.32 ± 0.83 


-5.77 ± 0.76 


3 


00 


± 





25 


(16) 


18000 


± 


900 (16) 


lo.stli 


150 (16) 


HIP 70275 


241 ± 6 (2) 


4.24 ± 0.69 


-10.74 ± 0.64 


3 


64 


± 





20 


(1) 


22790 


± 


1140 (1) 


11 3+^-1 


68 (2) 


BD -2 3766 


41 ± 10 (4) 


1.77 ± 0.59 


17.26 ± 0.62 


3 


70 


± 





20 


(4) 


22000 


± 


1000 (4) 


10 i + i-'' 


200 (16) 


PB 5418 


152 ± 10 (4) 


8.25 ± 1.48 


-4.74 ± 1.63 


4 


00 


± 





20 


(4) 


19310 


± 


966 (1) 




150 (10) 


HIP 56322 


254 ± 9 (17) 


3.54 ± 0.73 


11.52 ± 0.69 


3 


57 


± 





25 


(1) 


25501 


± 


1275 (1) 




160 (16) 


Ton S 308 


89 ± 10 (4) 


-0.38 ± 1.57 


-0.96 ± 1.76 


4 


00 


± 





20 


(4) 


17821 


± 


891 (1) 


6 0+"-* 
°-"-0.7 


120 (4) 


PHL 2018 


108 ± 10 (4) 


3.87 ± 1.19 


-11.27 ± 1.17 


3 


70 


± 





20 


(4) 


19095 


± 


955 (1) 


'■'-0.9 


150 (4) 


HIP 52906 


84 ± 11 (17) 


-7.39 ± 0.61 


1.34 ± 0.60 


3 


39 


± 





15 


(1) 


18898 


± 


945 (1) 


9 0+1-2 


160 (16) 
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Table B3. Computed values for stars belonging to group A. The values shown are: the present coordinates of the star {X, Y, Z) (in 
galactocentric Cartesian coordinates), the distance, the evolutionary age, the flight time, the ejection velocity and the coordinates (Xq, Vq) 
of the point of ejection in the galactic plane. 



Name 


X 


Y 


z 


d 




Age 


^flight 




^0 


^0 






(kpc) 


(kpc) 


(kpc) 


(kpc) 




(Myr) 


(Myr) 


{kms~^) 


(kpc) 


(kpc) 


EC 19337-6743 


-7 60+"-^" 


-0 -^o+o io 

"■^"-0.10 


-0 ■«i+°-io 

"-•^"-0.10 


0.56i2 


30 
.22 


r,r,o 0+314-4 
J32 


21-013:8 


28111 


-r. r,7 + 0.B3 
■'-■' ' -0.62 


-4, 


7,1 + 0.74 
— 1.05 


EC 05515-6231 


-8.00_Q 


,-1 ,-(1+0.20 
-0-60-0.30 


-0 40+°10 
"-^"-0.10 


0.71+0 


29 
.22 


36.7t«|;i 


31.0;-;-? 


3511° 


-4 78+"-°* 
'"-0.80 


-7 


-03+°-?* 
— 0.95 


EC 20089-5659 


-7 40+°-3'' 


-0 20+" °° 
"-^"-0.10 


-0 40+° l" 
"-*"-0.20 


0.81+0 


35 
.25 




25.9j:|;°o 


53 + 21 
^■3_14 


-4 21+"-77 
*-2'--0.74 


-4 


.90+°-=? 
— 0.67 


EC 03462-6813 


-8-UU_o 00 


-0 60+0-2° 
U.DU_o.20 


-0 60+°-l" 
"-""-0.20 


0.83+0 


27 
.21 


1 0+75.3 
l-**-1.8 


TO + 2 - 2 
1Z.^_2 J 


7q+20 
'"-16 


-6 59+"-07 
"-^"-0.50 


-3, 


-20+?-°3 
— 0.66 


EC 06682-5816 


o ,„-l-0.10 
-».1U_0 00 


-1 10+0-2° 
-'--'"-0.40 


-0 60+° l" 
"-""-0.30 


1.3lt"o 


43 
.33 


29 0+13-* 
^"-"-28.4 


35.81?-* 


1621=1 


-2-02l=2:=2^ 


-3, 


or) + 0.52 
— 0.39 


EC 10087-1411 


-».2U_o JO 


-0 90+°'^° 
"-»"-0.30 


o.6olS;fo 


i-ioi"o 


36 
.27 


09 .^+37. 8 
' -32.7 


18-71?:^ 


Mllf. 


-6.481J-J? 


-2, 


-84+°-™ 
— 0.64 


EC 06490-4510 


-».au_o JO 


-1 40+0-30 


-o-8otg:fo 


1.74^0 


.56 
.43 


18.3_jg-3 


23.711-1 




-6^47i;-.™ 


-6, 


.03+;-?^ 


EC 10649-2953 


-7 90+°-i° 

' ■»o_Q 00 


-1 iSO+0-30 


o.8oi°o:fo 


1.7ll°o 


.55 
.42 


36 4+*2-7 


42.9l?-| 


56llt 


-o-sii^-.S 


— 7, 


.86+?-S 

— 0.64 


EC 19071-7643 


-6 70+"-^° 
°' '"-0.30 


-1 20+0-30 
'-^"-0.40 


-0 90+0-20 
"■""-0.40 


2.03+0 


75 
.56 




6 6+°-^ 
"-"-0.3 


2261^3 


c ,1 +0.91 
^-21-0.68 


-2 


-,0+0.15 
■ —0.18 


EC 19476-4109 


n 0/1 + 0.70 
-6-30-0.50 


00+" °° 

"-""-o.io 


-0 90+0-20 
"■""-0.40 


1.93tS 


77 


37.lt^22:? 


34.3+^-5 


111+** 


-1 98+"-75 
-"°-0.82 


-4 


-6i+!;-s 

-0.64 


EC 20292-2414 


-6 40+°-^" 
°-*"-0.40 


60+"'^° 
U.DU_o.20 


-1 00+°-^" 
l-""-0.40 


1.97+;; 


'67 
.49 


7 4+*-6 
'-*-7.4 


16-211:^5 


13ltf6 


-5 92+"-=° 
"-"2-0.49 


-1 


cc+0.67 
- — 0.51 


EC 13139-1851 


'■^"-0.20 


-0 90+"-2° 
"-""-0.30 


1 10+0-40 
l-'"-0.30 


1.65i"o 


54 
.41 


1 7+24.5 
l*-~-18.7 


26.llt-« 


1221*3=2 


-4 14+1-51 
*-l^-1.28 


-4 


03+°-** 
-"•'-0.44 


EC 2()14()-(i935 


rn + 0.50 


-1 00+"-^° 
'-""-0.40 


-1 20 + 0-30 
l-'^"-0.40 


2.14iS 


79 


22 4+1-0 


12 + °-** 


1141*? 


-5 7K+0-76 
■'- '"-0.56 


-3 


qO + °-*= 
--'"-0.52 


EC 09414-1325 


^ „,, + 0.20 
-*^-'"-0.30 


-2 -^O+O-OO 
'^-■^"-0.70 


, .10 + 0.50 
l--^"-0.30 


2.78+S 


91 
.70 


1 + 6I.6 
1-^-1.2 


35.711+* 


105lfg 


-5.011?-,=? 


-6 


14+1-02 
-'^-0.94 


EC ()324()-(i229 


-7.8U_o 00 


-1 40+0-30 
'-*"-0.50 


_i t-o + 0.40 
l-°"-0.60 


2.12+2 


.73 
.54 


140.0t|3-^ 


45.51^:^4 


76111 


-0.98llj-7? 


-8 


oQ + 0.72 
■•'"-0.81 


EC 20153-6731 


-6 10+°™ 
°'^"-0.50 


-1 20+0-30 
'-^"-0.40 


>-,o+o-«'o 

l-""-0.50 




.96 
.68 


rjrj .,+ 17.9 

"- ' -15.4 


OQ 1 +2.8 
-^-^■^-2.4 


154l^° 


-1 52+l-°l 
^-"2_o 92 


-3 


7, + 0.66 
■ '2-0.59 


EC 19586-3823 


^■•^"-0.70 


10+0-10 
"-'^"-0.00 


-1 60 + °-'''" 
'--""-0.60 


3.201J- 


.25 
.91 


,„ 0+I6.8 
•'"■■'-II. 9 


1 +5.5 


1531^3 


1 93+"-83 
^■"•5-0. 75 


-6 


,4+0.73 
- ' -1.00 


EC 06387-8045 


-6 70+"-*'' 
°' '"-0.30 


_„ ifi+O.SO 

•^-i"-i.io 


-1 70+°-*" 
1- '"-0.60 


3.75tJ 


.29 
.94 


12 4+3-0 
'^^■^-9.0 


24.4l=-"o 


134+17 

Jg 


-3 -10+1-00 
•'-1"-1.47 


— 5 


,r+0.56 

■^■'-o.es 


PG 1533-1-467 


-7 70+0.10 
' ■ '"-0.00 


, ,0+0.20 
l.ou_o 20 


1 «n + °-30 
l-**0_o.30 


2.2-i+_o 


38 
.33 


C.+9.6 

"■■'-o.s 




213+37 
"^'■'-30 


-4 •<2+"-77 
^-•'2-0. 65 


-2 


.,0 + 0.19 
■■'"-0.17 


EC 20252-3137 


-5 40+1-"" 
°'*"-0.70 


c;n+0-20 

"-'^"-o.io 


1 sn+0-50 

-i.SU_o go 


3.2ltJ 


,24 
.86 


I2.5t«;| 


45.7+r-^ 


1551^3 


, 77+0.65 
~2- ' '-0.S4 


-4 


95+0 08 
■■^■'-0.80 


EC; 05438-4741 


-8.90_j, 20 


_„ i,, + 0.80 


-1 90+0-50 
'^■""-0.60 


3.67^0 


.26 
.91 




35.5113,3 


1001^8 


f. 0,1+1.79 
-6-82-1.19 


-8 


07+1-39 
■"'-1.62 


EC 20104-2944 


-4 70+1-3'' 
'"-1.00 


70+0-30 
"- '"-0.20 


-1 90+°'^" 
'--""-0.80 


3.90+1 


.59 
.14 


96 7+51.2 
' -26.7 


17-9li.^8 


172lf3 


-5 87+1-3* 
"-"' -0.99 


-3 


-1^+0.94 
-l''-1.29 


EC 06616-6107 


-8.00_o 00 


-3 40+0-**° 
•^-^"-1.20 


-2 00+°--''" 
^-""-0.70 


4.0lllj 


33 
.00 


4.2l°;22 


1 3 Q + 3-2 
io.y_2 3 


ISSlfg 


_4 51+1.84 
*-"l-1.17 


-6, 


02+1-2* 
"2-1.69 


PG 1205-H228 


jj ,„+0.00 
-8.20_o JO 


-0 30+0 00 
"-^"-0.20 


2 20 + °-^" 
^-^"-0.60 


2.22+0 


.76 
.56 


18.212+1 


14 4+*-^ 
-3.7 


254+0° 


-5 42+2-°l 
0-*2_j Jg 


-2, 


Og+0.28 
■""-0.17 


EC 23169-2235 


-7 30+0-20 
/.ou_Q 20 


n «n+0.20 
O.6OI0.2O 


, on + 0.60 
-^-•3U_o.70 


2.47tg 


.82 
.61 


+ 33-3 
"■"-0.0 


20.413-0 


15113*, 


-5 39+1-27 
"■•'"-0.95 


-4, 


on+0.86 

■80-1.02 


EC 09462-1403 


q en+0-40 


-4 00+1-00 
^■""-1.30 


2 40+0-SO 


4.89^1 


.61 
.20 


1 4+''2-l 


96 1 +13-5 


272102 


-4 86+3-*2 
't.oo_ J 90 


-3, 


,7+107 
' 1 22 


PG 2345-H241 


-9 00+°-^° 


■3-OU-0.90 


qn+0-70 

-^.au_j 00 


4.93+j 


.65 
.25 


, , ,+ 18.5 
ll-''-11.2 


26.7l5-t 


1491?? 


q ,,+0.79 
"■^'-1.30 


-3 


o„+0:93 
-80-1.14 


PHL 169 


-6 00+0 *0 
"■""-0.30 


„ qn+0-60 
•^-•^"-0.50 


-3 00+0-50 
■^■""-0.50 


4.87l°o 


.88 
.74 


3i^st^:*3 


21.4l|,0^ 


1881=20 


Q o,+0.63 
"-•'2_o.92 


-0, 


07+0.99 
■"' -0.82 


EC 10600-1358 


S 411+0-10 
-8.40_o 20 


-4 20+1-1° 
*-^"-1.50 


„ co+1-20 
■^■^"-0.90 


5.46+1 


.93 
.38 


68.7tlj=3;° 


55.1_j^3_3 


224l|3 


04+=-77 
"■"^-2.80 


-4, 


jg + 2.23 
■1°-1.21 


PG 1511-H367 


-6 90+0-20 
°--^"-o.io 


1 90+0-30 
±.yu_o 30 


3 60 + 0-™ 
■^■""-0.50 


4.27±"o 


72 
.63 


,„ c+20.8 
zo.o_22.7 


20.51^2:2 


256132 


-10 22+0-71 
lU..i.i_0.98 


0.37l°:?? 


EC 20011-5005 


-1 40+2-^0 
-'-■^"-l.go 


_, O0+0-40 


-4 20+1-^" 
^-^"-1.60 


7.87+-3 


12 
.24 


-16 2 + 31-1 
l°-^-16.2 


34.711732 


274137 


^■"'^-1.59 


-3.18ll-f, 


BD -15 115 


o on+O-OO 
-».JU_(j JO 


1 00+"-^° 
l-""-0.20 


-4 60+°-™ 
*-°"-0.90 


4.73iS 


87 
.72 


25 4+5-1 


30 4+^-7 
•'"■^-2.6 


249l*3^ 


1 7S+l-*2 
-1-78-1.24 


„ 77+0.68 
■'- ' ' -0.63 


PG 2219-1-094 


-6 30+0-30 
"'■■^O-0.30 


5 60+1-°° 

0-OU_o go 


-4 60+°-™ 
1-OO-0.90 




37 
.11 


25 3+*-'' 
^='-■'-4.0 


39 3+11^0 


■^""-40 


— 1 4 07+2-57 
l*-"'-3.97 


, ,0+2. 62 
^■2»_2.69 


HS 1914+7139 


-10 50+o-*o 


10 60+1-^° 
1U.DU_ J 50 


4 sn+0-80 

4.SU_o.7o 


ii.87j:ij-°'; 


„ 7+9.2 
2^- ' -11.7 


58 3 + 54.5 


,,,+ 175 
^1^-84 


-6 37+8-00 
"-^'-7.78 


-2 


gO + S.02 
■""-12.79 


PG 0855-1-294 


-13.80tl-j° 


-1 70+0-*0 
-'- '"-0.60 


4 90+1 **" 
^■""-1.20 


7.81+j 


80 
.99 


34.6tp^ 


24 1 +3-8 


278lf2 


-16 03+2-08 
'"-"■'-3.98 


-2 


gg+1.07 

■""-0.68 


PG 0955-1-291 


-11.60+;-?« 


-1 30+0-30 
-'-■^"-0.40 


4 90+1-*" 
^■""-1.30 


6.141; 


.87 
.51 


85.5+25;7 


7„ ,1 + 48.0 


2771*3! 


„ ,,7+11.26 
8-27_4 92 


-1 


07+2.58 
■"'-1.76 


PG 2229-f099 


-6 50+0-50 
°-''"-0.40 


6 00+^-°° 
o.uu_ J 50 


-5 10+1-3" 
^■-'"-1.70 


7.98+j 


.70 
.92 


34 5+11-0 
■5«-O_20.5 


63 l+*3-8 
"^-1-19.6 


1781S 


-3 35 + *-23 
■'-^"-5.05 


-7 


05 + 5-02 

-"•'-o.os 


PG 1610-1-239 


-4 00+°-™ 
"'-""-0.60 


3 40+0-00 
'^•*"-0.50 


5 20+°-^" 
^■^"-0.70 


7.45^; 


34 
.15 


62 1 + 13" 


^2-"-i3.i 


212lf3 


f,», + 2.50 


-3 


56+3-71 
-""-6.95 


PG 2111-1-023 


-2 40 + 0-»0 
^-^"-0.80 


7 50+1-20 
'-'^"-i.io 


-5 30+°-™ 
■'-•'"-0.90 


10.86ilj-^;5 


81.0 + lj0-3 


IKK K+**-2 
186.6„j 


24ll2| 


-6.58+3-?? 


— .5 


09 + 3-78 
-""-2.30 


PG 0009-1-036 


„ „„ + 0.20 
-8.90_o.30 


., r,o+1.10 
•^-""-0.90 


-5 60+1-*" 
■'-""-1.90 


6.6712 


16 
.65 


34-21;,^- 


40 s + 20.3 
^"-■'-10.7 


256 _ 40 


-5.91 + 3-72 
— 2 . b b 


-9 


q +4.16 
--'•'-8.64 


PG 0122-H214 


-12 50+0-™ 
lZ.OU_(j go 


4 70+0-90 
*- '"-0.70 


r c:n + 0.80 

-5.60_j 00 


8.62+J 


55 
.29 




39.9l^-^ 


156+32? 


-12.93llj-7? 


-3 


oQ + 2.02 
-•'"-2.14 


SB 367 


7 "in+o-io 

-/.»U_(j JO 


n sn+o io 
-O.8O_0.20 


-5 90+°-**" 

O.MU_j JO 


6. 02+ + 


.03 
.89 




ere, + 8.0 

oo.o_g 5 


198l*« 


-5 42+1-72 

"-*2_ J gg 


-6, 


g4+2.01 
■" -2.06 


EC 04420-1908 


-14 70+2-00 
'"-2.70 


-5 10+1-^0 
0.1U_2 JO 


-6 20+1-0" 
"-^"-2.50 


10 44+*-35 
-'"-^^-3.16 




29.3ll7g0 


257l°l 


-11 43+*-13 


-7 


40 + 2-74 
-*"-4.29 


PG 1332-H137 


-6 20+0-^0 


-0 60+"-^° 
"-""-0.20 


6.30+^j-™ 


6.54+2 


.13 
.70 


CO , + 13.4 

a».^_jj g 


32.4l?:| 


413+?? 


, or,+ 1.22 
"-■'"-0.98 


1 99+*-08 
j-.yy_2 22 


EC 19696-5356 


3 20+" °° 


-3 10+0-so 
■^■i"-i.io 


-7 30+1-00 
' ■■'"-2.60 


13 Rl+*-*° 
ld.»i_3 g3 


33.5_ jg 5 


29.2l^°jl 


4751^3 


-1 67+*-*l 
l-°'-4.01 


-6, 


00+4.10 
-08-6.77 


PHL 346 


-4 60+° ™ 


3 oo+o-^o 

■^■""-0.50 


-7 30+1-2° 
'-^"-1.40 


8.55 ^ 33 




28.5+|-7 


418lt? 


-3.891?-=! 


4-281?:?? 


PG 0914-HOOl 


— IS Qn+2-^0 


-13.8013-70 


10 90+3-00 
iu.yu_2 90 


20.621°-^* 


73.8^17 ? 


«.oi«i-i 


3691^5? 


-14.82lll7?3 


-24.72+?2=7^ 


PG 1209-H263 


-ii.8oi°o:fo 


-2 90+"-"° 
■'■""-0.80 


30.70t7;io 


30.93t;-,°J 


76 l+32'6 
'"--'-20.1 


192.6l57;2 


qqo + 293 

■'""-100 


— 26 7o+-'^5. 18 
M. 'U_gg g2 


-32 44+3°-87 
^-— 68.29 
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Table B4. Computed values for stars belonging to group B. The values shown are: the present coordinates of the star {X, Y, Z) (in 
galactocentric Cartesian coordinates), the distance, the evolutionary age, the flight time, the ejection velocity and the coordinates (Xq, Yb) 
of the point of ejection in the galax;tic plane. 



Name 


X 




Y 


z 


d 






Age 


*flight 


"ojec 


^0 


Ya 




(kpc) 




(kpc) 


(kpc) 


(kpc) 




(Myr) 


(Myr) 


(kms 1) 


(kpc) 


(kpc) 


HIP 81153 


-7.00+Ji 


20 


0.20+°=° 


0.50+°-™ 


1. 


1 0+1 
■I8I0 


20 


4.2 + 2-1 


6.8l|^3 


169+;2 


_ ss+0^39 

/ . 00 Q Qg 


-0.79+55-;* 


HIP 1241 


-8.101° 


.10 


0.30+g;° 


-0.60+°- 1° 


0. 


-64lo 


'16 


174-8+tJ'l 


29.112-J 


43+1 


-6.27+°-|° 


-6.47+°-|J 


HIP 28132 


-9.ooj:° 


!20 


-110+030 


-0.60+°-™ 


1, 


,59lo 


■37 


43-5+^5% 


19.9+3-J 


58+12 


-6.42+S'lo 


_ ,,_+0.S6 

0.00 Qg 


HIP 2702 
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